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SYNOPSIS 


COMPARATIVE STUDY OP STRUCTURAL BEHAVIOURS AND 
DESIGN OP LONG CYLINDRICAL CONCRETE SHELLS ' 
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for the degree of 
MASTER OF TECHNOLOGY 
By 

PRABIR CHANDRA BASU 
to the 
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Indian Institute of Technology, 
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The thesis presents the comparative study of the 
structural behaviour and design of long cylindrical concrete 
shell having span to radius ratio greater than 3 , axe assumed 
as long shells and such shells with edge beam, are considered 
for the investigation. Both reinforced and prestressed concrete 
shells are studied. Only edge beams of prestressed shells 
were provided with prestressing force. The response of the 
stress resultants and deformations are taken as structural 
behaviour parameters. 


Pgramatric study was made to investigate the influence 
of different shell parameters and prestressing on the behaviour 
evaluated by both the theories and to find out the domain of 
applicability of the beam theory* It was fbunid that inr case of 
reinforced concrete shell^beam theory is applicable for all 
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av^d 

normal values of radius* to thickness^ depth, of .edge beam to 
span ratios, semicentral angle and span to radii®-- ratio greater ; 
than 3. , within 4 to 15 percent of error. j 

Prestressing reduces the stress resultants and 
deformations with a suitable combination of depth of edge beam, 
and *rt44ec centric it ies of the tendon. The discrepency in the 
assumptions of beam theory increases with the increase of depth 
of edge beam and prestressing force. Restraints against the 
horizontal and rotational deformations yield. the stress resultant 
which are in nature to compensate each other. No shell with 
coefficient of prestressing more or equal' to 0.4 and depth .of 

edge beam to span ratio more equal to 0.04 could be designed by 

. . ' ' ' j 

beam theory. j 

that I 

It was found /the arch analysis by column analogy method; 

. ' I 

considering a thin strip of the shell with edge beam is not j 

efficient. A simplified formulation for arch analysis based on }; 

strain energy method, satisfying the compatibility of deformation; 

at the interface of the shell and edge beam, is developed. With 

the help of this arch analysis a systematic and fairly simple j * 

method, detailed in six steps, is presented. The proposed method; 

gives good results for single barrel shell, 

A series of table is given to make the computational 
work easy for the proposed method. 



CHAPTER-I 


INTRODUCTION 

1*1 GENERAL: 

"" " ,> 4 1 -: " 

Shells are one of the economical, effective and 
elegant structural system for the purpose of roofing of 
a large unobstructed area. They tend to carry loads 
primarily by direct stress acting_.in their plane, because 
of their geometry and small flexural rigidity of skin which 
cateigorise them to stressed skin structures. Among the 
different types of shell^oyllindrical shell is most ^popular 
for roofing purpose. They have extensive application in 
constructing airport hangers , industrial buildings, audito- 
riums, exhibition stalls, godowns, markets, garage units, 
bus terminus etc. 

For small span, Shells are.- not provided with edge 
beam, but as span increases edge beams are provided to 

reduce the stresses. For very long spah prestressing of edge 
beam yields the optimal utilization of structural properties 
of cylindrical shell roofs. Prestressing of edge beam has 
certain advantages as it reduces stress resultant and hence 
offers greater safety against buckling; deflections; 
c q asuBSfptio ns of steel and dead weight of the 'Structures* 

With the suitable selection of end and mid eccentricity and 
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magnitude of prestressing force it is sometimes possible 
to eleminate transverse moment altogether so that stresses 
in shell reduces to those corresponding to the membrane plus 
a constant,- compression. The cracks due to shrinkage, 
temperature variation can be tackled very efficiently by it. 

1 . 2 STRESS RESULTANTS : 

In general the stresses of a continunac 'sftr.uc.twral 
elements are shown in Pig. 1€ and they are , U jrty , 

^ * For a uni ~ ax i~ 

symmetric shell = N0 X . >■ )tor rejzzlindrical shell 

» 

Q x , » have very insignificant magnitude. N x , .. 

N 0 , and are considered as the main stress resultants in 
this investigation. 

1.3 LITERATURE REVIEW : 

Considerable work had been done on the analysis 
and design of Cylindrical Shell. However discussion here 
is restricted to the work available on analysis and design 
of ‘‘-cylindrical shell with prestressed edge beam by classical 
and 'Beam Theory'. 

In 194-9 H. lundgren( 3 )-> first Introduced a simple 
theory to analyse long .•,f$:lindrical shell by ’Beam Theory’. 

He compared the transverse section of a Cylindrical shell as 
a curved section of a beam and analyse the shell for N x and 
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N x 0 considering the shell as a beam spanning ir or. 0:10 cage 
tr a m to another in longitudinal direction,:’t'or I ‘'’tha 
calculation of and he performed the arch analysis 
considering a unit width strip of shell p s an arch loaded 
by the external loading and which is calculated by 

beam calculation. The &3eeh calculation that he did was 
for the interior shell and analysed it by column analogy method. 
J. Chinn (12) in 1959 demonstrated that the interior shell 
of a group of shells without edge beam even with l/R >1.11 

y 

can be analysed by beam theory, for urbh calculation he 
applied the elastic centre method. He also compared stress 
resultants obtained by beam theory and ASCE Manual (31). 

In the same year A.L ^.rme and H.W. Conner (14) introduced 
tables for the correction of stress resultants obtained ty 
beam theory. Ram&swami (7) described ■ hJmost all of the 
so called exact theory to analyse -cylindrical shell. He 
also demonstrated the procedure of designing the --cylindrical 
shell with prestressed edge beam and analysed amnterior 
shell of a group of shell without edge beam having L/R=2 

using column analogy method for arch calculation* G-ibson(l) 

in 196® presented the method to analyse the cyllindrical 

a 

shell by assuming it as/sha3p-ow she^l*. A detailed method to 
analyse a shell with pre stressed edge beam was put by him 

and a"' comparison of 3 J X obtained by exact method and beam 
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theory for shell without edge beam and with prestressed and 
non prestressed edge beam had been done by him. Tix'nual 
No. 31 of ASCE on the Design of Shell Roofs (27) contains a 
series of comprehensive tables ar-aa charts which reduces the 
lot of computational hazard. P. Dayaratnam(8) illustrated 
completely the procedure to design a single barrel shell 

withi 1 prestressed edge beam using ASCE Manual No. 31. 

D.P. Billington(2) discussed the classical shell theory by 
theory of shallow shell and use of ASCE Manual-31 and 
simplified the design problem. He also demonstrated the 
reason of non agreement of stress resultant obtained by exact 
solution and beam theory for a single barrel shell without 
edge beam of L/R-2.5 and introduced a chart for stress 

resultant prepared by be am. arch approximation. Marshall(9) 
had shown it is sometimes possible to eliminate transverse 

moment altogether so that the stresses in the shell area reduced 
to those corresponding a membrane plus a constant compression. 
Ramaswami and Saeed(34) used Beam Method to analyse a 
T^pa.rb'dlb id ' purb'o 1 odd shell . 

Very little work is available on' the structural 
behaviour of the cylindrical shells with edge beams by b$«l» _ 
theory £.or reinforced or prestressed shell. 
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1,4 PROPOSED IWESTIGATIQF: 

Design of cylindrical shell by so called classical 
method is complicated and involves tedious computations, 
which^as a result ^js n ot an attractive procedure to the 
engineers.-in practice. 3eam theory,- sufficiently accurate 
for reinforced concfete cylindrical shell, involves less 
complicacy. Though the method described in ASCE Manual-31 
is not that complicated as other classical methods ' hut ‘if is 
not suitable for a shell with any type of combination of 
L/R , R/t and $ K . Again prestressing is desirable for 
optimal utilization of stressed skin structures. So there 
is a need for research on beam theory to explore its • 
applicability to analylyso the stress resultant of long 
cylindrical shell with prestressed edge beam. 

The summary of proposed investigation is: 

1. Analysis of long cylindrical shell with edge beam 

by classical and beam theories. Then make a comparative 
study of the structural behaviour. 

2. Parametric study of the structural behaviour of the 
shell by 'classical and beam theories. 

3. Develop the formulations for ar-ch. analysis by strain 
energy method and illustrate the procedure. 

4. Series of tables are given to design barrel shell roof 
by the above method. 
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1 .5 CO-ORDINATE SYSTEMS AHILSI fi tt CQWEl€lQM S ; 

— M l— I H I , m mm rnl l mm n Sm mmmm i — 

The co-ordinate system adopted in the investigation 
is shown in the Figure 1.1(b). The following sign conventions 
are -used unless otherwise mentioned (ref. fig. 2.2). 



N x = longitudinal force considered positive if it is tensile. 

= Transverse inplaneforce in outward' - direction is 
positive. 

■®xj# = Tangential shearing force is positive when is • creates 
tension in the increasing direction of ft and x. 

M.0 = Transverse moment is positive if it induces tension at 

the bottom fiber. 

u x = longitudinal deformation is positive if it increases 
along positive x-direction. 
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tiy = It is positive if it increases along positive 
y-direction« 

u r = Radial displacement is positive if it acts outwrxds. 
w v = Vertical deformation in downward direction is positive. 

= Horizontal deformation in out war'd! direction is positive. 
w x = Longitudinal deformation is positive when it is 
exte*sional' in nature. 

w m = Rotational deformation is positive when it has anticlock 
wise sense. 

For other- purpose general sign convention will be 
followed. Throughout the study sign convention,- as mentioned 
above, will be followed unless otheE-wise mentioned. 

1.6 CLASSICAL SHELL THEORY ( 8 ) : 

Classical shell theory developedyhas the following 
steps or system, 

t 

i) Brimary system:- Shell resisting surface load 
solely by direct inplane stress resultant i.e.- membranec state. 

i 

ii) Errors:- Forces generated by membrane' - : theory at , 
edges which does not satisfy. tho' - bounisu"y conditions. 

iii) Corrections:- Corrective line loads axe applied 
p.t free edges to sat isf^ -ictu/' Lboundary conditions. 

iv) Sending Solution:- These corrective line loads 
will cause bending of the shell for which bending solution is 
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required. Membrane- solution is obtained considering, that, the 
shell is loaded at the edge by corrective line loads. 

.Complete solution yields by superimposing membranes solution, 
obtained considering shell is loaded by ext ernal loads, on. ben- 
ding solution, considering edge disturbance as the loading 
on the shell. 

a) Assumptions : 

1. Materials are homogeneous, isotropic and obey 
loolce *s~ lav; . 

2. Deformation of shell produced due to loads are 

small. 

3. Thickness of shell is small when compared with is 
is other dimensions and radii of curvature. 

s v 

4. Kirchoff's hypothesis holds good i.e., 

i) Rectilinear elements normal to the middle 
surface of the shell remain rectilinear and normal even after 
deformation and do not change their lengths. 

ii) formal stresses acting on the planes parallel 
to the middle surface are neglected inf c Q ifitn turner 
stresses. 

5. Poisson* s ratio of concrete is neglected. 

1 

6i Shell is assumed to be simply supported 
between two longitudinal trover cos. 
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b) Membrane Theory : 

Membrane theory assumes external loading and the- 
surface of this shell are continuous. 

The equilibrium equations for membrane state of stress 
(ref. Fig. 1.&;)) are as follows: 


m 


x 


1 




3x 


R 


9)2 


+ q R * 0 
x 


R 


9^ 

~~w 


3E 


xj z5 


+ qyR = 0 


9x 


( 1.1 ) 


( 1 . 2 ) 


+ q/ =0 


(1.3) 


Membrane solution is obtained by integrating above solution 
which are given below, 

i) For loading acting on the chord area of the shell 

(q u ) 

12 2 o , . 5 nroc 


N * R ( — ) ( Cos 2 0* - Sin 2 0*) —y 

x R xr% 


Sin 


1 

(1*4) 


1 3 nirx , 

R-jjC = - q R ( - 5 -) Cos 0' Sin 0' Sdos (1.5) 

F u R nn 1 


^ = - q u R Cos 2 0'±S±n 


nmx 


( 1 . 6 ) 



ii) For the loading action on the surface of the 


shell (q g ) 


1 s 2 


R x = - q g R (-£-) Cos 


^7 


“ < 3 * ( 4 "> Sin0 ‘ 


nTc 


NjZf = - q g R Cos jtf'^Sin 

/ 


mtx 


nnx 

Sin — j 

(V 7) 

n.7tx 

(1.8) 

■ Cos 


(1.9) 


The expression of membrane displacements are given below, 
i) For q u , 

R 1 3 p 

' ' 2 a ft 


u 


In ^ 


x Et ' R 


( )'(Cos^ 0* - Sin 2 0* ) - T .» Cos 

R nr 


mtx 


1 

(1.10) 


u = - q R JL. ( J— ) 4 Cos 0* Sin 0 f 2+ ( 

7 u Et R L 


R 1 4 


mtR . n7tx 

) / Sin 


"r = - 1u R g£— (, ff— > 


» nrcR 4- 2 ^,1 6 

■v ( ~ — > 003 0 [ w Sin 





(1.11) 

Sin 2 

0' ) 

o< 

rntR ; 

-- L ~ j 

i 

6 

Sin 

nrcx 

(1.12) 
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ii) Eor q , 

s 


u = q g R 


R 


Et ' R 


u 


u. 


R 




R 

Et 


( ) ^Cos 0* Sin 


nux 


3T3 


(1.13) 


■n it 


1 4* o o R 2 

^ On/-> ^TT. f Y^-in (X * 


- 1« R S" ( — : (^-)“)Sin *' j-f Sin 


nTtx 


. 1 

(1.14) 


4 4 

, 1 A 2 2 R. 2 n. 7i ft ,4. 

( ) (1+2n %< ) + — ( ) ) 

R k 1 * 2 1 


ZCOS0< 


Sin 


mcx 


(1.15) 


n 1 %' 


where , 0' = 0 X - 0 

Of. loading 

Eor the expression 1. 4- to 1.15 the sinusoidal variation ^along 
the longitudinal direction is assumed. 

° ^ Sending Solution : 

The equilibrium equations for an element on the shell 
(Eig. 1.2), Are as follows,' 

9N-, 


’x 9N x0 


9x 


-it + 


d0 


+ R = 0 


(1.16) 


9% 9N x rt 

-j-j — + — R-Q0 + 10R-O 


(1.17) 


X 
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9Q *_ r + SQ 0 


3x 


a*5 

"a ~ 


3x 


+ E0 + Q. R = 0 


3M. 


X0 R + Qtf R = 0 

a x F 


R - 


9M 0: 


+ Q R =0 

30 x 


(1.18) 

(1*19) 

( 1 . 20 ) 


In the above sets of equations total number' of .equations 
is 6 and number of unknowns is 9. To solve this statically 

indeterminate systems the other conditions are taken from 
stress-displacement relationship which are given below; 


U x = Et 

% 

^x0 


Et ( 


9u x 

3x 

1 


3u t 


u. 


R 


w 


R 


Et 

2~ 


3u- 


7 


3u-> 


3x 


30 


( 1 . 21 ) 

( 1 . 22 ) 

(1.23) 


M. 


Ef 


2 

a u r 


1 I 


x 


12 


Et' 


ax* 


3u y 


Et 


12 

3 


R* 


M x0 = 


12 


2R 


30 

9l3 y 

3x 


1 

"4* - 

R^ 

1 


2 

a u r 

"a0 2 ~ 


) 


3x30 


(1.24) 


(1.25) 


(1.26V 
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with, the help of the equations 1.16 to 1.26 the final 8th 
order differential equation with respect to w r obtained 
as y 


8 12 a^Up 

v W "ta* 5 


12 rVi Zd3%fS 


1 


St 3 


1 

T 


.2 


Where ^ 


9x90^' ' 

1 3 


R 


0 *27) 


r — 


ax" -bT 


\ . « , 


Homogeneous solution of eqn. (1.27) has the form 


u r = S' 


UTC 


n=1 ,3 


k^e Sin kx, where k = -y 


(vt^8) 


putting Eq. (1.28) on the homogeneous' part of 1.27 the 


characteristic^ eqn. stands as follows; 
.2 , 12 


r 2 ( - h 2 + —g ) 4 + 


k 4 = o 


(1.29) 


x Roots of the Eq. (1.29) 


m y^_ = i ( a + ip ) 


** ± ( a* + ip ) 


(1.30) 

(1.31) 


where 


\ 


n ( i+"”y ) ' v + i + ( 1+ Tf ) 


(1.32) 
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hi ( 1 + v y + 1 - ( i+ t ) 


a ' 

■P*' 

V = 

S’ 

V 


\ 


{ ( i-V ) 2+ i - ( i-TT ) 


f ( i-V ) 2 +i +< i-- f ) 


(iL ) 2 
\ 


R 2 


3 (*R ) ( _) 


(1.33) 


(1.34) 

(1.35) 


(1.36 ) 

(1.36) 


Therefore , 

u r =pL 1 e (<, 1 +iP 1 ) * + A 2 e ( “l * 1|3 1 > * 

- + A 5 e ( “l' + «i*> * + A 4 e (a i - * 

+ + * + kJ- * (1 ’ 37) 

5 6 

+ A^" a 1 + 1|3 i^ ^ + A s e^" “l " 1 ^l ^ J Sin tat 


Expressing A-j Ag in the form r 0-+ il) and 

measuring the angle 0 ' from the crown the stress resultant and 


the displacement can be expressed by following expression, 

E = x(a B 1 - bJ3 2 ) ( a + bB 1 ) | 2 + ( cB^ _ dB 4 ) f 3 

- ( eB 4 + 4B 5 ) (1.38) 
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where f 


= Cos (3 ft* Cosh a ft* 
Tf 2 = Sin p ft* Sinh a ft* 

^ = Cos p * ft* Cosh a’ ft* 


( 1 . 39 ) 


Sin |3 ’ Sinh a*\ fti 


for different stress resultant and displacements, the shell 
coeff ib-iehts- hre give.npiii fcafcle A # 1 # 1 and A.1 .2 in Appendix-1 

Particular solution:- The particular integral of Eq, 1.27 
may be taken as 


q = Cop ( - 0 ) Sin kx 

X 


(1.4-0) 


Substituting of the Eq. 1.40 in Eq. 1.27 'gives 


C = 


Where q 


2 A 2 .2 

R k + 4 k + 2/R 

( R^t 2 /12 ) ( k 2 + 1 /R 2 )V 4 


(1.41 ) 


4 

q g 

71 & 


If the loading. is q u then the particular solution is 


<l x = - % k*’ where 1 = q u 

tt; 
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1.7 BEAM THEORY : 

Analysis of long cylindrical shell by fees.® theory 
has essentially two steps, i) beam analysis , and ii) arch 
analysis. 


a) Assumptions : 

Beam' theory is based on the same assumptions as 
classical theory given in artiole 1.6(a), except that 
certain factors are assumed to be negligible and they are: 

i) Relative displacements within each transverse 
cross sections. 

ii) Longitudinal bending moments M^. and radial shear 

forces Qx on the shell. 

iii) Torsional moments ( M x ^ ) 

iv) Strain from inplane shearing force. 


b) Beam Analysis 


1 


The object of beam analysis is to calculate B x , 

O.G. of the section of the section is: 

Crawn-rf-g A c = 2 R$ k t (1.42) 

Distance of the C*G. of the section 
from the crown is ; 

* i 1- Sin 0 k / j^) R 



V* 1 

ffigk 1 .3 ' Curved Beam Section 


(1 « 43 ) 
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Moment of "irgrtia of the section is; 

I c = k I R3t » k l = + Sin 0^ (Cos 0 k _ 2 Sin 0^/0^ ) 

(1.44) 

First moment of area at any section 0 is; 

2 0 

q = k Q R -t , k Q = ( Sin 0 - Sin ^ ) (1.45) 

4 

Distance of the top fiber from C.G-. = Z^_ = Z c 
Distance of the bottom fiber from C.G. = = h - Z Q 

Therefore 


R x at top = q 0 l 2 Zy ^' 

(Comp. ) 


N x at bottom - q Q 3 / 
(tension) 

R x 0 section 


t Z. 


'yb 


lx X 

T T (1 ’ “T" } 

ix x - 

rr (1 - — > 


( 1 . 46 ) 

(1.47) 


q. >1 
H o 


Q 


l x 

{- ) ( 1 . 48 ) 

' 2 1 7 


c) Act ch Analysis ; jU- \)Hs- 

Object of the aycfii analysis is to calculate , M^c 


and 


An elementary ai^cfe- strip, which is under the action 
of !u» % and specific shear, is considered. Eq. ( 1 . 45)express 

th.e magnitude of at a section 0 along x. 



Therefore , 


(1.49) 


OT x0 _ VQ 

3x 2 I 

So the specific shear at any longitudinal section is constant 
and ‘- its vertical component will act upwards and horizontal 
component inwards. The ihiiLQh is divided into a no. of nodes of 



Vertical projection of ds = dZ^ = R ( Cos ^+1 “Cos./#^) 

( 1 . 50 ) 

Horizontal projection of ds= dy^ = R(Sin$}-Sin0^ + .j ) 

(1.51) 

Therefore, = q u dy* + q ds-^ + dzj, (1.52) 

6 3x 
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Therefore , 
i 

M $i = .2 


3=1,2 

i 


^ Y .i) + F hi^ z i ** 2^) 


hj v i 


(1.54) 


~ T #i - ( 1 D . ) Sin 0’ - ( E ' F h1 ) Cos 0? 

3=1,2 J 1 3=1,2 ' J 1 


(1.55) 


■'i ' 3=1,2 

1 ,8 DISCUSSION: 


( L- 3rd > Cos ~ ( J =1>2 Sin (l *56) 


The beundary conditions for arch calculations 
is determined from the physical system of the shell itself, 
which can be divided mainly in 3 types which are given below; 

1) Single Barrel:- The arch may be considered as axisymertr ic 
statically determinate system, loaded by external loads and 
supported by specific shear, having two ends -freer , 



Dig. 1.5 ARCH MODELS 

2) Interior Shell:- The arch may be considered as axisym'ni&tric 
statically inderminate system, loaded^as, type-1 , with both 
ends fixed. 

3) End Shell:- The arch may be considered as axi-asymetric 
statically determinate system, loaded as type 1, with one 
end fixed and another end free. 

In the next chapter both the methods are illustrated 


by numerical examples* 



CkAPIER-XI 


ANALYSIS OP GYLIIIDRICAL SHELL BY CLASSICAL AMD BEAM THEORIES 

2.1 imODUOTIO U ; 

The stress resultants and the displacements of a shell 
can be expressed through four constants a , b t c , d which can 
be determined from boundary conditions. For the shell without 
edge bean the boundary conditions are that , llg 

are zero at the free edge-. If the shell is provided with edge 
beam, the boundary conditions are the edge deformations of shell 
are same as those of edge beam at the interface of the shell 
and beam. 

2.2 THEORY O F ED Q-E BEAM; 

For the shell with edge beam the compatibility of 
deformations are: 

1 . Compatibility of vertical deformation 

2. ' Compatibility of horizontal deformation 

3. Rotational compatibility and 

4. longitudinal deformation compatibility. 

If the stress resultants and the deformations are denoted by 
the subscript ’b' then, 

Vertical deformation, w.y- g =* w rkSin 

Horizontal deformation w^ s = Cos ~+Pyk (2.1) 
Rotational deformation w ms = 0^ 

Longitudinal deformation w xs .»* 




22 


Loading of the edge beam at the interface is the edge 
stress resultants of the shell edgd in same magnitude but of 
opposite sign along with the external loads, if any, and self 
weight of the edge beam. If V , H , II and S are the vertical, 
horizontal, rotational and shear loads respectively at the edge 
of the shell, the edge deformations are given by: 

w -vb = w vb < v > H > M > s > lb- ) ' 

w hb = w hb ( V , H , M , S , q b ) 

%b = "mb ( V . H , M , S , q b ) (2 ' 2) 

w xb ~ w xb ( V , H , M , S , q b ) 

^here q b represents the external loading on the edge beam. 

Again V , H , M , S , w vs , w^g , w mg and w xs are 
associated with a , b , c , d discussed in previous chapter 
and they can be expressed as functions of a , b , c , d if the 
loading and geometrical properties of shell and edge beam are 
known. Hence eqs. of 2.1 and 2.2 yield the interface compatibi- 
lity as follows , 

W vs ( a, b , c , d , q s ) = w yb ( a , b , c , d , q b ) 

Whs ( a , b , c , d , <ls) = w hb C a » b » c * d » % ) (2.3) 

w ms ( a , b , c , d , q s )= w mb ( a , b , c , d , q b ) 

w xs ( a , b , c, d, q s ) = w xb ( a , b , c , d, q b ) 

m ■' ■ 7 7 : ■ ■■ f: | 
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where q b and q g are the total loading system on shell arid 
edge beam. The four equations of 2.3 will now determine 
the four constants. 


2.3 DEFORMATIONS QP EDG-E ffEAM: 

Details of loading on edge beam are shown in Pig. 2*1. 
The expression of deformations are given below; 



w vb = ’V'T + w vH + w vM + w v 

Where the first subscript indicates the direction' of 
deflection and the second indicates the force causing it. 

b) Longitudinal Deformation? 
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Wxii r - JT $7^7 e M = =b33»* (2.6) 

w xW = 0 

w xb = w 3*. + w xS + w xH 

<5 ) Horizontal Deformation: 

*! 15“ 2 „ TT 

w hY = ~ . ~2~ — e h e v 7 = e b31' y 

2G b J 

VhS = ■ it; Sh s = 01,528 

WUH = — 2“ I -4“ -1-+ ” " 2 e T e hl H = %35 H (2*7) 

E b I z 2G b J. 

1 % *— 2 

w hM = “ — 3r e v M = Cv,2r M 

m 2 2G b J v b 34 


w hb = w hV + w hS + w )iH. + w hM 
d) Rotational Deformation: 


w nrV = 


1 

~"7~ 

1 

TT- e v Y 

> 

& 

O 

11 





w mS = 

0 






1 

V-2 

" °b43 H 

w mII “ 


m 2 

— ’ e v 2 

2»b J 

%H = 

1 

— 

n - 2 M 

2 G b J 

°b44 M 

w mb *“ 

w m¥ 

+ WmH + w mM 





Prestressed cable in a vertical plane will cause 


deformation in vertical plane only. 


'These deforu. 


nations ■ are vertical and honiqitud&aideformations. The edge 
beams are pre stressed by curved cable (tendons ) haying parabolic 
profile. The effect of prestressing the edge beams is 
equivalent to the particular cases of loading which are 
described below and shown in Pig, 2,2, 

1. Continuous distributed upward pressure( w^g) 
due to the change -of curvature of the cable profile and if 
magnitude is given by, 
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Case 1 


M 0 (Case 3) 



Case 2 and 5?', 


Pig. 2.2 EFFECT OP PRESTRESSI¥G FORCE 


2. End moment due to eccentricity of anchorage is 


given by, 
\ = 


^ 0 0 


( 2 . 10 ) 


3. Axial thrust applied at both ends ®f the centre 

4 

of gravity of the beam section.given by , axial thrust = — . 

% 

( 2 . 11 ) 

Vertical deformation caused by prestressing force, 
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i eg l 2 - 

m4 n ^' 


e o , Pe 
®bfy l2 


( 2 . 12 ) 


and longitudinal deformation, 


w 


xp 


_4_ 

it 


( b-+ 

^b 


e o®v 


E b^y 


m 


8g e v 

"TT“ 


E b Jy 


•)p e (2.13) 


The coefficient 4 -/it is used in all the Eqs. 2.12 and 2.13 as a 
Fourier coefficient of uniformly distributed load converted 
into sine series. 


2.3 ANALYSIS OF SHELL BY CLASSICAL THEORY: 

The theory described in article 1.6 is the theoretical 
aspect of cylindrical shell analysis by the method described 
in ASCE Manual-31. A series of comprehensive tables are given in 
ASCS Manual-31 for barrel roof analysis. In this investigation 
shells are analysed using those tables taking n=1 . 

All the coefficients given in the table of ASCE Manual-31 for 
the loadS) assuming sinusoidal variation of loads along the 
longitudinal direction. The method to analyse the shell is 
given below, 

Membr ane Analy sis:- 

= Ar 1 RC 1m Sin m (2 ' U) 

= -i- ( g— -J q R o 2m C*9 mx 
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Nj# = <1 R c 3m Sin mx 

Where q is the loading. For q = q g , c 1m , , c^ m are 

given in column 7 , 8 , 9 of Table 1 B and for q = q^ column 1 , 

2 , 3 gives the above coefficients. If c g = lV(R-^t-E s ), 

then the membrane displacement can be obtained from the following 

expression, 


i ) For q 


wym = q g R c 


w hm ^g ^ c s 


. € 

/ 2R \2 — / \*r 

( __ ) + — „ + ( y c 


n 1 


— I 

7t^ 


R. \4 

1 ' ^wlm 


R 4 

(— ) c w2m Sin m x 


Sin mx 
(2.15) 


'wlm 


and c W £m are given in column 10 and 11 of Table 1 B. 


ii) For q = q u 


w. 


v.m 


^uR c s 




R /I T +( n R/l ) 4 /l2 




' wlm 


Sin mx 


w- 


hm ~ ^u®- 0 


( fc/l) 4 c w2m1 + <>1+ (icR/ 1) 2 /2 + (7tR/l)Vl2^X. 


c w2m2j Sin mx (2.16) 

c v;1m , c w 2 mi and c W 2 m 2 are given in column 4 , 5 and 6 of 
Table 1 B and the longitudinal deformation, 

= —L- ( Nx) / (1/R) 2 = f osm (2*17) 

^ E s t E s t 

Errors and Correction:- To satisfy the boundary condition 

i.e. E x , , Ngf a i shell edge are zero, a set of line loads 
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should be applied at the edge which is equal in magnitude 
but opposite in sign of and N >^k* Correction 

forces are resolved into vertical and horizontal components 
and they are : 


V L = - %ic Sin ft Sin m x 
H l = - Cos 0 Sin m x 

S L = - Cos m x 

B ending Analys i s : - 


<2.18) 


The correction line log,ds will disturb the membrane 
state of shell and bending of shell will occur. The bending 
analysis of shell is done in following way. 


Edge deformations due to correction line loads 

Sin mx 

H t jsin mx. (2.19) 


3 

& 

1} 

o 

— — i 

o 

<J 

-A 

+ <V-2 S L + c v3 

1 2 1 r 


Wxl " ( R ) tE s [ 

c n v i + c i 22 

r 2 r 


Wml " 5^ L Cmi 

Y L + c m2 s l + C 1 

w h 1 = C h 1 V 1 + 

°h2 S 1 + c h3 H L 


c]si: 


m3 


’13 

Hj, J Sin mx 


Where c v1 , c v5 , c^ 2 .. c hl , c h5 , c n1 , c m3 and c m2 are given 
in column 1 y fL, & and 9 of Table 2B and c^/e^ 

c -| 2 are given in column 1 , 5 and 9 of Table 2JL. 

When the shell is provided with edge beam a set of ' 
forces, say v^, H b ,.s t and will generate at the inter- 
face. The deformations due to the interfabe forces are: 
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w vbi 

w xbl 

%bl 


r 


G s jjVl 7 b + c v2 s b + c v3 H b + c v4 


] 


Sin mx 


= c 


Ta.bi = 


(_1 /R )2 1/Est [_ C H Y b+- ®i2 s b + c i5 H b + c u M b /a 
s j^ c h.1 v b + c h2 s b + G h3 H b + G h4 M bA ] Sin mx 


R 2 

E T 
s A 


C m1 7 b + c m2 ^b + c m3 *% + c 


3 m4 


M b /R~ 

.. 


Sin mx 
( 2 . 20 ) 
Sin mx 


where c m1 , c^., c m2 , c v4 , c h4 and c m4 are given in column 
3 , 6 , 9 , 10 , 11 and 12 of Table 2 B and c u in column 4 
in Table 2 A. If positive interface load acts on the sKell 

th§i negetive load will act on the beam, then the deformation 
compatibility, at mid span of shell edge yields, 


w 

m 

+ w vl 

+ w, , - 

vt>i 

w vb 

= 0 


w 

xm 

+ w xL 

+ w xb\ - 

w xb ■ 

= 0 

(2.21) 

w hm 

+ WhL 

+ w hbl ' 

w hb 

= 0 



w ml 

+ w mbl “ 

w mb 

= 0 



( 3rd Subscript 1 indicates deformations otuo.- to JoccdU/Y^) 

c s c v4 


c s c v1 + c b1t °s c v2 + G b12 c 


s c v3 + c b13 


(1/R) 2 J. Ci-t+c, 


R 


E s 


■^11 +c b21 


✓In 2 1 /In 2 1 " ,1% 2 c 14l 

U> c 12 +c b22 — c i 3 +c b 23 %> 

s 


R Egt 


E s t 


c s C h1 + c b3i c s c h2 +c b32 c s c h3 + c b33 


; s c h4 


ix 


_ / 
B- R' 

v" m1 ° M1 v 


G m2 


R‘ 


E„I 

S x 


C m3 +C M3 


R 


c s c m4 


,+c 


b34 


R 


«+ o 


b44 


P 


b 


•tel 


& 


(2.221) ! 


2 

l°3 


!C 


4 
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Where , 

°1 = w vm - ’Vt + w vbg + w vbp 
°2 = -w sm - w xt + w xbg + w xbp 
C 3 = ~ w hm ~ ™ht and c 4 = - «ml ' (2.23) 

w vbg == vertical deformation due to self weight of the beam 

= c b11 b d f . (2.24) 

w xbg = longitudinal deformation due to the self weight of the 
beam = c b 2i b ^ T (2.25) 


Analysis of shell by using the formulae direct from the 
solution of differential equation of shell involves the 
evaluation of the constants a b , c and d as basic unknowns 
whereas method using the tables takes the interface forces 


as basic knowns, which are determined from the boundery 
conditions. Once V b > s b > * and M b are known ^the total 


correction loads ^ 

V = V L + s b , S = ^ + s b , 

Stress resultants due to bending; 


H = H L + H b , M = % + M b 
(2.26) 


= ( * 1t V + C 12 , S + c 13 H + c 14 ) x ( l/R) 2 

R 

R x 0= (c 21 Y + c 22 S + c 23 H + c 24 JL_ ) x (x/ R ) 


S 0 = c 51 Y + c 32 S + c 33 H ■+ c 34 M/a (2.27) 

= ( c 41 Y + c 42 S + c 43 H + c 44 M/R)R 
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Total stress resultants can be obtained by super -imposing 
membrane stress resultants obtained from Eqs. of 2.14- to 
bending stress resultant obtained from Eqs. of 2,27. 

It is general practice for single barrel shell to 
assume that edge of the shell does not offer any restraint against 
horizontal and rotational deformations, so Eqns. of 2.22 
reduces to 2 simultaneous equations. For the interior 
shell of a multi barrel group * 0 and w m = 0 while for 
single shell H = 0 M = 0 . Numerical computations of the 
method is illustrated by the fallowing example, 

2.6 EXAMPLE 1 

1 = 40 m , R = 12 m , t = 0.12m , h - 0.15 m , d = 2.0 m_ 

= 40° > e Q 0.5 m , c m = 0.7 m , p e = 400,000 kg, 

loading : 

p 

i) Live load on shell = 75 kg/m (acting on surface of 

shell ) 

ii) Load from lime terracing = 45 kg/m 2 (acting on surface 

of shell) 

iii) Self weight of shell - 2400 x 0,12 = 36 0 kg/m 2 

Total load on the shell = 480 kg/m 2 

Grade of concrete of shell and beam is M350, 

OVng. SiW&U is 


t 
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a) Membrane Analy si s : - 

Membrane stresses are calculated from Eqs. of 2.14 
and' are given in Table 2.1. 

Table 2.1 MEMBRANE STRESS RESULTANTS 


0 

(kg/m) 

(kg/m) 

(kg/m ) 

(kg/m ) 

40 

- 1 4040 

o . 

-6234 

0 

30 

-13820 

-2297 

-6139 

0 

20 

-13190 

-4524 

-5858 

0 

10 

-12160 

-6614 

-5399 

0 

0 

-10750 

-8503 

-4775 

0 


b ) Errors and Corrections : 

Errors are = - 4775 Eg/m and N = - 8503 kg/m . 
The vertical and horizontal component shear correction line 
loads are calculated from Eqs.. of 2.14. 

= 3070 kg/m 

H l = 3658 kg/m 
S L = 3503 kg/m 

c) Bending Analysis: 

i) Deformation of shell edge at mid span:- It is 
general practice to ignore the restraint' offered by t;he edge 
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beam towards the horizontal and rotational deformation as the 
beam is slender one. So for this problem H b .and' H b may 
be taken as zero. From Eqs. 2.15 , 2.17 and 2.19. 


W_ + Wi T -t 

vm vi 


w xm + W x1 


c s 54300 


1 o 1 

( r 26130 

* V 


ii) Deformation of Edge Beam:- The contribution of e^ 

towards the beam deformation is very small and hence ignored* 
From Eqs. pf 2.. 5 , 

b w vbs = C b1 ( ~ 0,12319 ) V b + — j— 0.19351 S b 


From Eqs. of 2.25 

w vbg = C b1 112 * 3 


From Eq. 2.12 

w Tb p - e b1 1.2543 e m - 0. 2923 e 0 
From Eqn of 2.6 


w xvb + w xbs = . c b2 ^ °* 60793 ) 7 b - _ 1.2732 S b 

w xbg . . = c b2 557 * 3 
Erom Eq. 2.13 


w xbp 
Where , 




4 


c b1 : 


bE b d : 


and c b £ 


bE b d* 


iii) Compatibility of Deformations:- Vertical deformation 
compatibility from Eq. 2.21 , 
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1 R ,3 <1 

42.33 V b + 0. 9561 S b - — ) (-0.1231 \ + — 0.19351 S b 

+ 112.9 + ( 1.2548 e m - 0.2923 e 0 ) p e A 2 ) + 54500 = 0 

(2.28) 

Longitudinal compatibility from Eq. 2,21 

-0.216 Y b - 0.1668 S b - Ji- ( JL- ) 2 (0.60793 V b - JL_0. 12732 S b 

D CL 1 U 

Pe 

- 557.3 + ( 1.2732 d t 1.4478 e 0 - 6.1923 e m ) ) ■+ 26130 = 0 

L (2.29) 

d ) Total Stress Resultants ; 

Eqs. 2.28 and 2.29 gives V b * - 1226 and S b - - 15930 
Therefore, final correction line loads from cqs. of 2.26, 

\ V = 1844 , S = - 7477 and H = 3658 

Bending stress resultants are calculated from Eqs. of 2.27 &&&■ 
adding these with membrane solution total stress resultants 

are obtained and given in Table 2.2. 

Table 2.2 TOTAL STRESS RESULT ARTS 



*x 

( kg/m ) 

• 

(kg/m ) 

*0 ' ' 
(kg/m) 

M 0 

(kgm/m) .... 


(at x=l/2 ) 

(at x = 0 ) 

( at x= 1/2) 

(at x=l/2) 

40 

-51 180 

0 

-5393 

735 

30 

-43180 

-10100 

, -4919 

695 

20 

-22950 

-17090 

-3674 

991 

10 

-27290 

-19610 

-2100 

865 

0 

-44570 . 

-1.9720 

-619 

0 
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Longitudinal stress at the bottom of the beam 


If. 


xb 


1 * 

b d 2 


0.60793 ( i-b d f - \ ) + 0.6336 ~ S b 


% 


L. 

4 

+ (_d + 1 .4478 e 0 - 6.1923 e* ) 


( 2 . 30 ) 


297000 ( kg/m 2 ) 


The variation of stress resultants are shown in Fig.’ 2,5 (a) 
and (b). 

2.7 ANALYSIS OF SHELL BY BEAM THEORY: 


Theoretical aspect of beam theory has been described 
in the article! .7. The beam method may be devided into two parts 
as beam and arch analysis. To illustrate the numerical compuia- 
tations of the method the same shell of example t • is analysed 
below. 

a) Beam Analysis : 

The transverse section of the shell along with the 
edge beam is assumed to iconstiute an mncracked 

section of a beam whose geometrical properties are shown in 
Fig. 2.3 » 

Chord length, = 2R Sin = 15.4286 m 

Depth of the section, h = R ( 1- Cos ) + d = 4.8075 m 

Area of the Section, A s = 2 ( Aq + Ap ) = 2.6106 m 2 
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Distance of the C.G. from crown, Z c = 


Moment of inertia of the section, I g ; 


•^c ^rc ^ 


2(A c Z c +A b Z b ) 


1.6077 m 


2 ( I b + I 0 + A b Z^ b + 
5.3950 m + 



Fig. 2.3 GEOMETRICAL PROPERTIES OF ASSUMED BEAM SECTION OF 
EXAMPLE 1 


First moment of area of any section, 

i) for shell portion, Q = 2Rt Z s - R (1-Sin fajfa ) 

. • n 

p 2 

ii) for beam portion, Q = / 2b Zdz = b(Zy b - Z ) 

Co-ordinates of any point on the shell, 

Z t = R ( 1- Cos ) - Z 8 

Y ± = R Sin fa 

Total load U.D.L. on the shell on the chord area, \ 

i 0 = f U s y + 9g ) fV sln + 4al X 0 = 8277.4946 kg/m 
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Bending moment ;at the mid span, 

% l2 3 

M = = 1655.2 x 1CK kg. m. 

8 

Shear force at, the end section, 

Itl . , 

S. = 165.52 x 10 ? kg. 

1 2 3 



The section is a symmetric section, so the stresses of the 
half section is given in the Table 2.5 
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Sable 2.3 HT X AND OP EX. 2 (BY BEAM METHOD ) 


Hodal 

Points 


Co-ordinates 

Q 

(kg/m2) 

(at x= D/2) 

(kg/m) 

(at x=B/2 ) 

Z. 

i 

Y i 

1 


• 3.2028 

7.7136 

0.000 

107700 

0000 

2 

- 

2.2028 

7.7136 

0.8099 


5218 

3 

- 

1.2028 

7.7136 

1*320 


8504 

4 

40 

1.2028 

7.7136 

1.320 


8504 

5 

35 

0.5644 

6.8832 

1.540 


9921 

6 

30 

0.0028 

6.0 . . 

1.6093 


10370 

7 

25 

-0.4804 

5.0712 

1.5461 


9963 

8 

20 

-0.8816 

4.1040 

1.3732 


8845 

9 

15 

-1.196 

3.1056 

1,1089 


7146 

10 

10 

-1.4228 

2.0832 

0.7775 


. 5010 

1 1 

5 

-1.5596 

1.0464 

0.4002 


2570 

12 

0 

-1.6077 

.0.0000 

0.0000 

702635 

0000 





\ ■■ 




b ) Arch Analysis : 

.fin elementry arch strip of width <tx i s considered 

■ ■ ' . • : ' ; ■■ * 1 . ■ . ■ ■ ' : ' I ... 

at midspan, whiph containts the shell and -the edge beam 

together and divided into 12 nodal points as shown in Pig-are 2. 

The details of arch calculations are given below in tabular 

far*m Table 2.4 contains ds. , dy , da , __^x^ # Table 2.5 

9 x 




40 


contains 


Qx S 

and Table 2.6 contains N^ 


8N X 0 9N X $ 

dz , ds , q u dy , F v , 5^ (= dy ) 


0X 


, distance of the nodes 
from the top. It is assumed that nodes 4 and 12 have half 

1 

of the segmental length and 1 c", 3 have segmental length d 

4 

where 2 has d/2. 


Table 2.4 ARCH ANALYSIS 


Nodes 

♦ 

3N xi * 

3x 

ds 

ay 

dz 

1 

. ' . §, ■ 

000.0000 

0.5 000 

0,00000 

0.50000 

2 

260.8898 

1 . 0000 

0.00000 

1 .00000 

3 

425.7712 

0.5000 

0.00000 

0.50000 

4 

424.7712 

0.5236 

0.88230 

0.39528 

5 

494.9012 

1.0472 

0.85752 

0.60048 

6 

517.4841 

1.0472 

0.90660 

0.52344 

7 

526.5362 

1.0472 

0.94884 

0.44244 

8 

441.0083 

1 . 0472 

0.98364 

0.35808 

9 

356.7043 

1.0472 

1.01120 

0.27096 

10 

250.0121 

1.0472 

1.03920 

0.18168 

11 

128.9848 

1 .0472 

1.04280 

0.09132 . , 

12 

000.0000 

0.5236 

0.52334 

0.01140 
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Table 2.5' ARCH ANALYSIS 


* — 1 
o 

pj 

CD 

m 

4 a s 

(D 

w u dy 
( 2 ) 

3x 

( 3 ) 

( 1 )+( 2 ^+( 3 ) 

' ^ 

( Wrf Uy ) 

3x 

i 

213.6288 

0 

ooo.ooo© 

213.6288 

000.0000 

2 

427.2576 

0 

- 11.7789 

415.4787 

000.0000 

3 

427.2576 

0 

- 45.4222 

381.8325 

000.0000 

4 

427.2576 

0 

- 96.6527 

330.6049 

207.0933 

5 

427.2576 

0 ■ 

- 157.9163 

269.3413 

424.3885 

6 

427.2576 

0 

- 232.9607 

194.2969 

469.1511 

7 

427.2576 

0 

- 270.8705 

156.3871 

499.5986 

8 

427.2576 

0 

- 297.1789 

130.0787 

434.7934 

9 

213.6288 

0 

- 167.64932 

74.7334 

360.6997 

10 

180.0000 

0 

- 212.8856 

- 32.8856 

257.7425 

11 

360.0000 

0 

- 260.8898 

99*1102 

134.5054 

12 - 

180.0000 

0 

000.0000 


000.0000 


w~ for shell 408 kg/fc 2 and for beam 360 kg/m 2 
s 

Wu=^ k • 
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Table 2.6 

AND M0 ( 

BY BEAM METHOD ) 

Nodes 


ii Ehl? 

M0- 
kg m/m 

1/2) (at x= 1/2 ) 

1 

4.8110 

0000 

0000.0000 

2 

3.8110 

0000 

0000,0000 

3 

2.8110 

0000 

0000.0000 

4 

2.8110 

-1430 

0000.0000 

5 

2.1726 

-1580 

1635.1483 

6 

1.5954 

-1804 

2961.8528 

7 

1.1248 

-2136 

4094.0868 

8 

0.7236 

-2387 

4997,2154 

9 

0,3560 

-2591 

5683.6779 

10 

0,1754 

-2712 

6 160. 6666 

11 

0.0456 

-2794 

6428.3211 

12 

0.0000 

-2887 

6508.4988 


2.8 DISCUSSION : 

It is evident from the articles 2,6 and 2.7 that there 
is no agreement between the stress resultants evaluated by beam 
method and classical method. A comparative study for the 
structural behaviour determined by both the theories is 
necessary to get the idea about the error involvement in 
beam theory. The influence of the different shell parameters, 
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prestressing force and eccentricities of the cable should be 
investigated. In next chapter(3) a comparative study of 
the structural behaviour of the shell of example 2 and 
parametric study are dona. 




Cro ion 


C) r* l M** by k>«<m 

theory 


a) fk <£ Nx* toy classical 
theory 




$ 0 $ - 12 * 10 ' 


d) N# <C M# 


by Classical 


t fU A ttOto Of STR.ES& R.E&uUTAtotS CS*^ 






CHAPTER- III 


COMPARATIVE STUDY OE STRUCTURAL BEHAVIOUR BY BEAM AMD CLASSICAL 

THEORIES 

3 .1 IMP ODUC T IQ H ; 

Eron figure 2 ^ 5 it is found that the variation of stress 

resultants of the example-1 determined in section 2 .6 by classical 

theory do not agree with the stress resultants determined by beam 

theory in section 2.7. The object of this chapter is to invest iga 

our 

the structural behavi-/ estimated by both the theories. 
Investigation by parana trie study is done to find out the influenc 
of different shell parameters and the prestressing force on the 
structrual behavior. 

3.2 A HA LYSIS OE SHELL IN EUIIY MEMBRANE STATE BY BEAM THEORY^ 2 \ i 

The shell is in fully membrane state when it offers 
resistance against the external forces by it's inplane stresses 
only. A cylindrical shell with circular generatrix attains fully 
membrane condition when it's transverse cross section is contim«o 
i«e.» the shape is circular ring and loading is continuous -with 
respect to its transverse cross section- 

The shell, to analyse by beam theory-, is assumed as a 
beam s panning inftferc direction of it f s axis of rotation with a 
hollow circular transverse cross section. 




4 * 



Fig. 3.1 GIRCHuSTrS^ CROSS ’ 

SECTION 

Shear force'S- = - J_ 2 % R a Cos mx 
x m g - 

Z , = R Cos ( * /2 - $) = - R Cos £f 


Loading on the shell = 2 n Hq f . Sin me 

3 '° 

Moment of Inertia, I = itE t 

V 

• First moment of area at any section 
jZf, = 2R 2 t Sin 0 

Bending moment at any lognitudinal sect! 
section, M = 2 n R q . „,1 sin mx 


Longitudinal stress, <jt 


2 X Rqg Sin mx 


Shear stress. 


?: = - 2^ Rq c 


_1_ Cos mx 


-R Cos $ 

•rt r ' 5 t * 

2 R 2 t Sin 0 

2 itt R 3 t 


(3.1 


Trans verse force, = ” dg Sin • (~ E Cos ^jj-) 

Substituting and = cT t and m_ = n% /l in 

eqn. of 3.1* 

BL = q~ R ( -Jr-) 2 Cos 0. — 2_ gin rLTrx 

8 H nV -I- 

5 ^=-q g E(— ) Sm*( Cos __ (3.2) 

= q g R Cos 0 Sin n % x 

If 0^ = n is substituted then eqns. 1 .7, 1 .8, 1 .9 would give exact 
same expressions like eqa. of 3*2 




Scale far .Mx«f 


i2Xio* -4x10* o -tx» ( 
seal* 


VAR.1ATI0M OF57RESS RESy/lTA UTS 

( £*1 3 Pe=o) 


Beam. 

Classical 
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3*3 COMPAteiyE STUDY OP THE STRTJOTTTRATi BMIAVIOm 

our 

By the term structural behavi-/ it is meant the response 
of the shell by stress resultants and the deformations. As the 
investigation mainly aims to explore the feasibility of the 
aplicability of beam theory to design cylindrical shell the 
following parameters are considered for structural behaviour: 
r x , ^ 0 , at iaid span and at support . 

a) Stress Resultants 

i ) \T V The variation of 53" obtaired by classical method is 

not linear while beam method gives linear distribution. Eor both 

the cases magnitude at the bottom of edge beam is greater than 

that of crown and the ratio of 0" at crown to that at bottom of 

x 

edge beam far classical method is 0,1765 and for beam method is 
0.1222. Again XT a t crown by exact method is 3*287 tima$greater 
than that by beam method and at the bottom of edge beam the ratio 
in 2 .805 . For the case of classical method magnitude decreases upto 

V 

the section in between 20°to 10° then bifurcate towards the 
increasing mode and attains it's maximum value at the bottom of 
edge beam but for the beam method it follows through out the sane 
path t Of attains the maximmat the bottom of edge beam. 


ii) JT^ The mode of the variation of along the depth of 
the edge beam is approximately similar with a difference that the 
rate of charge of IT x ^ at t he region where it attains it‘s maximum 
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of magnitude for the classical theory is smaller than that for 
the beam theory but the main difference is that the magnitude 
obtained by the beam method is 50 to 60$ less than-, that obtained 
by classical method. 

iii) : The mode of the variation of ^exhibits it*s tendency 
to decrease ..towards, the edge of shell both for classical and 
beam method-. At the crown obtained by the beam theory is 

0.1 945 times lesser than that obtained by classical theory while 
at the shell edge it is 2.318 times greater. 

iv ) : Sign of M^, obtained by both the theories is positive. 
The mode of variation shows the decrease in magnitude as it 
proceeds from crown and ultimately 0 at the shell edge. For 
classical theory magnitude decreases upto a section wlWre 0 is 
30° and then increases upto section 0 = 40°. From this point the 
magnitude again decreases and cone s to 0 at the bottom of the 
edge beam. Beam theory always gives the lows? value of than 
classical theory. 

b) Deformations 

The defrom t ions may be calculated by the classical 
theory from the eqns. of 2.19 by putting V,S,H,M in place of 
Vg* Sg» H 2 and where eqn. of 2.26 gives the values of V,§,H 
and M. The def ormations of the example 2 is calculated as 
follows, 
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w 

v 


■ *s «v1 T+ ' %2 s + 

— 2.4 ~ 

-5 j 42 , 38 z 1744 - 22,99 x 3568 + 0.9856 (-7477) 

BT tE L 
s 


= - 48 m.n. 




h 


J s c h1 V + e h2 S + G h3 


H 


w. 


-= — I 22.99 x 1744 - 13.59 x 3658 + 0,2931 (7477)1 

R^tE L J 

- 35 m.n 
2 


111 


R 


E I 


°n1 T + V s + V 1 


E 1 O .1245 x 1744 - 0.07595 x 3568 + 0.001413 (-7477)1 


V 


■0.017 radians. 


Deformations by beam theory are calculated as follow: 


w 

V 


i4 r 


it El 


,4 


7c4 


( ^t - 


5 P n 2 

L_» ) - 1 p a 

,2 %2 e e 


3 .5 mm 


'v^ and can not be determined by beam theory as it neglect those 
quantities by virtue of first assumptions of article 1.7(a). 

Beam theory gives the magnitude of w less than that of 
classical theory, as in this case w y obtained by classical theory 
is 13*21 times greater than that determined by beam theory in 
opposite direction. Beam theory gives and negligible 
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where as 


6 

vh 


for classical theory they are not negligible: actually 

= 0,728 and 5 = 1 ,9 

vm 


Where, 



(V 1 ) stained by classical theory 
(W^/l) obtained by classical theory 


6 

vm 


( ^ir/^k ) 0 bi a i*ied by classical theory 
(W y /l) obtained by classical theory 


3.4 lEFIOEHCE 0E HESTRESSIEG- FORCE 

4 

To investigate the influence^ presstr ess irg force , 

shell without prestressing force is considered as refer eic e« The 

stress distribution for the same shell of example -1. of article 

2.6 is shown in fig. 3.1 . It is found that the ^agreement of 

and obtained by beam and shell theory is close enough but 

for E^ and it differs. The contribution of P g and external 

loads along with the dead load of shell towards Q - ' is shown 

from fig* 3*2 it is observed that contribution of P , calculated 

0 

by exact theory, is more prominent near the region of shell edge 
& edge beam and less significant towards crown* while it's 


* ani it's variation is highly nonlinear along the depth of the 
transverse section of the shell-.* 
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influence follows linear path far beam theory. Both the method 
give-, contribution of P g as tensile stress at crown and 
compressive at bottom. The influence of P , evaluated by beam 
method, at the crown is 7.89 times greater than that, calculated 
by exact method while for the bottom of edge beam the ratio 
is 1 .81 . 

3 .5 EFFECT OP BOUTORI COIOITIOMS 

The boundary conditions and Uhe deformations of the 
shell, with and without prestressing force, determined both by 
beam and classical theory are given in the table 3.1 

Table 3.1 DEF0RMA.1 10 KS OF SHE1X ED&E 


Bound ery 

Assumption of 


Deformations 



condition 
for class- 
ical theory 

beam theory 

P =400,000 kg 

CD 

ii 

O 



Classical 

Beam Classical 

Beam 

w = W -v 

vs vb 

w,„ = w , 
vs vb 

(assumption 1 

of section 

1 .7(a)) 

= -48 mm 

3 ,5m.m 

21 mm 

15*6 

H = 0 

(w^ exists) 

w^=neglifeible 
(assumption 2 ) 

w-i = -35 m.m. n egli~ 
h ^ gible 

1 mm 

negli- 

gible 

M = 0 

(w m exists) 

w m = negligible 

w m~"’0*077 
• (rad ) 

negli- 

gible 

0 

negli- 

gible 
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For P g - 0, 8^^ _ 0.0476 ard = 0 . so for the shell 

without prestressing force w h and w m may be assumed a s negligible 

but with prestressing force they: can not be assumed as negligible 

as it is observed in the previous article that for P =400,000 

1 6 

^vh ~ °*728 aid 6 vm = 1 .9. Again for 

u /-> ( w v ) determined by classical theory 

x e = °» *7 — : — — S L Z LtQ - 1 i'345 

rained by beam theory 15*6 

while for no pries tr as sing it is 13.21. Therefore the assumptions 
for boundary conditions in classical and beam theory agre^when 
P Q = 0 which is not true far P g = 400,000 kg. 

3.6 PAR AM ICR IC STUDY 

our 

The structural behavt-/ of a cylindrical shell is very 
sensitive to the effect of certain parameters. Therefore it is 
desirable to study the influence of different parameters such as 
geometry of the shell, magnitude of prestressing force and 
eccentricities of the cables on the s tructural behaviours,. In this 
study the effect of parameters considered are, 

1 . Span to radius ratio, L/r 

2. Semi central angle, 0^. 

3 . Depth of edge beam, d 

4. Prestressing force, P 

G 

5. Ecentricities of the cable, e 0 and e m 
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Variation of one parameter is affected at a time while others are 
kept constants to observe it's influence on the stress resultants 
and the deformations* 


As the arch analysis is dependent on the beam analysis 
and if beam analysis yields accurate encugh longitudinal and shear 
stress resultants, then and could be determined approximately 
correct by proper arch analysis satisfying the boundary conditions. 
Therefore paramatric study is done for CT, and defromations .■ 

In the figures the variation of M^ is shown where longitudinal 
stress calculated by beam theory is close enough with that 
determined by classical theory. 

Tile non-dimensionalised parameters which are considered 
for the study are given in table 3*2 

Table 3 .2 BOB DIME NS 10 BA LIS ED PARAMETERS 


Shell paxaneters Mon Dimensionalised form Expression 


Semi-central angle 

4 

4 

Radious, R 

Span to radious ratio 

vk 

Depth of edge 

Depth to span ratio 

d/X 

beam, d 


T) 

Erestressing force 

Coefficient of prestressing p = 

Jr 

e 

o 

P e 

force, p 

ins. 1 

Eccentricities 

End eccentricity to depth ratio 

0 D 

O 


Mid eccentricity to depth ratio 

e n/ d 


<* 1 (T q u + q g V Sin 4 
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^ x * longitudinal stress, at the crown, decreases as 0 k 

increases from 30° to 45° in all I/R ratios but increases as l/R 

increases fur all values of 0 k for both p=0.4 and p=0„ At the 

our 

bottom of the edge beam the behavis*/ is just reverse. Both beam 

our 

and classical theory yields the above behavt**/ beam theory gives 
the value of at crown lower than that of classical theory 'for 
all values of L/R and 0 k when p=0.4 but for p=0 it always gives 
higher value than classical theory. At the bottom of the edge 
beam, bear, theory allways underestimates the value of iJH for P=0.4 
but for p=0 it gives higher value than classical theory for l/R=5 
and 0 k - 30, L/fe = 10 and 0 k = 30° and 40°. Rest of the cases it 
gives lower value than the classical theory. The agreement in the 
variation along the depth of the transverse section does not achiev 
for p = 0.4. When p = 0 the agreement is there and beam theory 
gives the varation of accurate enough except the two cases 

of L/R =3-33 when 0-, =40° and 45°. 0~ attains it’s minimum 

J£ X 

value in between the shell edge and 0 = 10° for all the cases 
except for 0 k = 30° and 40° when l/R = 10 and maximum value at 
the bottom of the edge beam by classical theory when p = 0.4. 

Por the same value of p beam theory gives minimum value at 

v 

the crown and maximum value at the bottom. (p x is of 
compressive nature by both the theory for p = 0.4 all along the 
depth. But for p = 0 maximum compressive stress occurs at 
the crown and maximum tension occurs at the 
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bottom of the edge beam by both beam ard classical theory. 

: ^ or rj ' 0n Distressing stress there is good agreement 
between beam and classical theory in the distribution of ISf^ , 

The minimum error is 3*5% while the magnitude of maximum error 
is 11%. For prestressing shell no comprasion exists in between 
beam and classical theories* For reinf creed shell error increases 
when smaller values of has been taken far higher Value of Tj/r . 

b) Deformations 

i) Vertical deformations: When p = 0*4 the vertical deformation:’ 
calculated by classical theory, gives the magnitudes much higher 
than those determined by beam theory. By beam theory difference 
between the vertical deformations far different values of is no, 
significant where for classical theory it is opposite. For p = 0 
mode of the variation is more or less same and the difference of 
magnitude is not much. 

ii) Horizontal deformations: Both for p = 0 and p = 04 can 
not be determined by beam theory while classical theory gives 
significant magnitude for both the cases. 

iii) Rotational deformations: Mode of variation far p = 0 is 
similar with that of w^. Ihrt when p = 0.4 it is to some extent 
diff erent»Classical theory gives significant magnitudes far both 
p = 0 and p = 0,4 while beam theory neglects it's magnitudes. 
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The variation of w^, and with l/R and are 
shown 'in. .'fig. 3,10. 

3-8 EFFECT OF THE DEPTH Of EDGE BEAM 

To study the effect of the edge beam on the structural 
behaviour, the values of d/L ratio consideration are 0 * 05 , 0 . 04 , 

0 * 03 , 0 . 02 . Other parameters are kept constant, as given, 

L/R = 0.03, $ k = 30°, E/t = 100, b/t = 1.25, P = 0,4 and 0, 
e o/d = 6 j/d = 0,4 

a) Stress Resultants 

The behavior of the stress resultants are shown in 

figure 3-5 and 3-6 aid the following observe! ions are made; 

i) KT • When p = 0 the agreement between the beam and the 
IK! ^ 

shell theory in the variation* is very much dominant for all the 

values of d/l. When d = 0.4 compression occurs both at the top 

and bottom by classical and beam theory where gr at crown has 

lower value than at the bottom of the edge beam for a value of 

d/L 0.5 to 0.4 classical theory gives tension at the edge of the 

shell at d/l = 0,5 and it’s magnitude decreases as d/l decreases 

and the shell is totally in compressive state when d/l =0.3. 

Beam theory does not give any where tension, - : - For d/l =0.02 

is 

the value of at the top/always less than at the bottom and 

X 

shell is fully in compressive state by both the theories, - The 
agreement in the variation/slightly achieved as d/l decreases. 
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: ®"° a S re e E ent between "he am and classical theories 

was found for prestressed shell# for reinforced concrete shell 
beam theory gives accurate enough stresses which is well comparable 
with classical theory. 

b) Def or nations 

Variation of deformations with d/L are given in fig #3 .10. 
from the study of this curve following observations are made. 

i) Vertical deformation (w y ); for p = 0 at d/L = 0.02 beam theory 

gives 0.722 time vertical deformation determined by classical 

theory but at d/L = 0.05 beam theory yields 1 .49 tines greater 

defroim tion than the classical theory. 7/hen p = 0.4, deformation 
is 

by beam theory/always less than the that determined by classical 

theory. As d/L increases deffereice between the defromation is 

negative 

increased. Classical theory always gives f deformation but for 
beam theory it is ^ negative for the range d/L = 0.04 to 0.05 . 

W y decreases for both p = 0 and p = 0.4 and magnitude far 'p = 0 
is always greater than the magnitude at p = 0.04 by both the theories. 

ii) Horizontal deformation (w fe ) : By beam theory horizontal 

deformation is negligible and it's magnitude can not be . determined . 

is 

w, determined by classical theory, for p = 0.04/always negative 
n 

and the magnitude increases as d/L increases, for p = 0 w^ is 
possitive and decreases d/L increases. At d/L = 0.0325 it is 0 
and ,f isr.. other values it is negative. 
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iii) Rotational deformation (w r ) t The variation of w with d/l 

is of similar type as that of w , . 

h 

3 .9 EFPECT OR ECCEMERTC ITY : 

The following combination of eccentricities are considered 
to study the effect of eccent icities in structural behavior. 

1 . eo/d = 0 2. oq/I = 0 3. eQ/d = 0.4 

V 4 = 0 V d = °- 4 e „/ a = 0 

4. ©q/I = 0,4 5. e o^ = 6. e 0^ = 

e n/ 4 = 0 * 4 . V' d = 0 - 4 ^ = 0 - 2 - 

other parameters are kept constants as follows 

Vfc = 3.33, d/I = 0.03, R/t = 100, b/t = 1 .25, p = 0.4 

our 

Prom thestudy of the behavi-/ curve for the stress resultants ’■ 
following observations are made. 

i) (f ; Por the combination 3 the agreement between beam and 

classical method is close. Por combination 1 defference between 

the CT determined by beam theory and classical theory is not 

much. Combination 6 exibits more-or less same manner of variations 

for combinations 2,4 and 5 no .agreement exists. Both the theories 

keep the shell in *the compressive state through out the transverse 

depth. Beam theory gives lower values of $~ x at the crown and 

bottom of edge beam than the values determined by classical theory. 

But CT" determined inthe region of mid depth by classical theory . 
IX 

has lower values than those determined by beam theory. The 
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influence of eccentricities on 
while at the. bottom of the 
inf luenee. 


l J at the crown have little 

•A. 

edge beam they have significant 


V * ^ or a -^ o ombinations, except 1, the difference 
between determined by the classical theory and the:: beam 
theory decreases as § E /d increases and Sq/i decreases# both the 
theories gives lower values of R when e^/d = 0»4 and for the 
values of §q/1 greater than the © m /d> both the theories gives 
higher values of Eccentricities have very sensetive influence 

on the magnitude of for both the theories. 


5. 10 EFFECT O f P RESTRESSIHj FORCE : 

Values of p considered for this study are 0.8, 0.6 , 
0.4, 0.2 and and other parameters are considered as follow: 
Ii/R=3 .33 &k = 30° , R/t - 10.0, b/t = 1 .25, d/l = 0.03. 
a ) Str es s Resultants: 

The behaviour curve for stress resultants are 
shown in Figure 3.8- From the study following observations 
are made. 

i) (T : Tor p = 0.8 beam theory give (T* at the 

x x 

crown as tension 'while for classical theory it is compresion 
for p = 0.2 and 0 both the theories give compression at the 
top and tension at the bottom. For other three values of p 
sr„ at bottom is compression and the magnitude decreasing with 
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when p is decreasing. Prom both the theory it is observed 
(J^ a "t crown increases as p decreases. The agreement between 
beam theory and classical theory is close for p = 0.2 and 0 . 

ii) 

The differences between the magnitude obtained by 
beam and classical theory decreases as p decreases and close 
agreement achieved for p = 0.2 and 0. 

b ) P ef or 1 ma tio ns : 

The behaviour curve shown in figure 3-11 . 

i) w v 

The variation of with p is linear for both 
the case. Beam theory gives the magnitude always less 
except when p = 0.2. for p = 0.4 w y has very small magnitude 
in posetive direction by beam theory while by classical theory 
its • magnitude is about 8 times greater in opposite direction. 

ii) v- 

The variation is linear. At p = 0 it is magnitude 
is not 0 but very smalland attains maximum magnitude at 
p = 0.8 by classical theory while beam theory assume it as 


negligible . 
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iii) w m :~ 

Same type of variation as that of only difference is it 
it is magnitude at p=0 is 0. 

DlSCIJSBIOSf : 

From the above investigations following observations 

are made, 

i) The discrepancy between the boundary conditions and 
the assumptions of beam theory is more significant for prestressed 
shell than non prestressed shell. 

ii) Beam theory can be applicable for add values 
of L/R , d/l and when p = 0 . 

iii) Beam theory may be applicable for a limited range 
of shell parameters for prestressed shell, 

iv) Restraints, fully or partially, yields the stress 
resultants whicfc compensates each other for 

In the next chapter on the basis of the observations, 
made from the above investigations, theory for arch analysis 
developed using casting liano's second theorem taking 
consideration the variation o± loads with Sin 
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CHAPTER -IV 


DESIGN OF CYLINDRICAL SHELL BY SIMPLIFIED BEAM A HD .ARCH AM LYSIS 
4,1 INTRODUCTION ; 

It is evident from Pigs. 3 .3 , 3-4> 3-5 and 3-7 that beam 
theory yields very good approximate stress distributions for 
N^ and N^ for almost all combinations of h/R , 0 ^ and d/l for 
non prestressed shell. But in article 3.4 it was found that the 
distributions for 1$ and 11$ by beam theory for single barrel shell 
do not agree with olassical theory • Again it is a well established 
fact that beam theory yields accurate enough stress distributions 
for N 0 Ifv , N 0 and 11$ for the interior shell of multibarrel shell 
group with and without edge beam when p=0. The object of this 
chapter is'- to establish such a formulation by which any type 
of cylindrical shell can be analysed and designed by beam and arch 
analysis , 

Like conventional beam method this method has two distinct 
parts, beam analysis and arch analysis, The beam analysis does 
not change but the arch analysis. is different. The theory of 
proposed method is developed assuming the loading on the shell 
varies with Sin n itx/l. 

In addition to the assumptions' made in article 1 -7, 
the method assumes, "The principle of superposition holds, 
in which case the strain energy is a homogenious function 
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of second degree, where the strain energy can be expressed as a 
function of external forces” * Eh is assumption is mainly for arch 
analysis » 

4 »2 BEAk A hA I if S IS 


The beam analysis is all most similar to that which was 
described in article 1 .7 (b). Ehe loadirg is assumed to vary along 
the longitudinal direction as half sine series. 


Loading on shell q. = Oq Sin ek 
W here 4 0 = 1„ ( 4 U + 4 q ) 


Si¥ s V 

Moment at ary longitudinal section, M = — r 

m* 


Sin mx 

Shear force at ary cection, S f = _1_ ^ 0 . 0 . ( ^ aad 

moment of inertia of, the sections are determined as before, only 
Qj_at any node i determined as follow, 

Qj_ = 2 E^t (Z s -R k zl ) 


At any longitudinal section, QH 


J yt 


1 


X 


At ary longitudinal section, 


s 

Q. 


q Q Sin inx 
(4.1) 


s 


m 

1 

-jp q 0 Cos mx 
(4.2) 


Q i 


3 x 


2 

q Q Sin mx (4.3) 


Specific shear, 


2 
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4-3 AMIESIS OF CIRCULAR jiBCB REST IBS OQLASglG POU1DATIOI 



Pig. 4.1 CIRCULAR ARCH OH ELASTIC FOUNDATION 


Moment at any section 0, 

M = VR (Sin0 k - 3in$ ) + EH (Cos^ - Cos#) + M s + ^ 

N = -V Sin# - H Cos# + 3SL 

S-f = -V Cos# - H Sin# + Sfg 

Where , K = Moment due to external loading 


3^ = Elastic redundant moment at support 

M - Axial force due to external load 

s 


Total strain energy of the arch. 


U « 


f 

0 


if 


2 31 


ds 


S 

+ /' 

0 


2 

W 


2AE 


ds 


(4-4) 
(4-5) 
(4*6 ) 
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From Qastiglianp ' s , second theorem. 


Vertical deformation at the support, w v1 


= a u 


6 V 


Horizontal deformation at the support = 3U 

3 H 

'Rotational deformation at the support, = fu 


(4.7 a) 
(4.7 b) 
(4.7 c) 




Bitting eqns. 4.4 and 4*5 on eqs» of 4.6, and equating w y1 , 

W m1 t0 w v> w h 81111 V 

+ Ik OT ) 7 + (T 3 !^) H + U?k mm > H* - - * 3 Kn (4 -S) 

< l3, W 7 - < l3k »ih + H + C^W *' = ®”h - I 3 ^1«K 22 (4.9) 

7 + 1* 3 ^) H + - 13 ‘Si (4 - 10) 


Where, 

k mv 


12^ k ( 2 Sin 2 ^ + 1 ) - 

6 Sin 2#^. 

(4 .1 1 a) 

^mh 

as 

12j^ k ( 2 Cos 2 ^ k + 1 ) - 

18 Sin 2j#k 

(4.11 b) 

^Snm 

3B 

CM 


(4.11 c) 

^xnvh 

«s 

12jZf k Sin 2jZ? k * 24 Sin% 


(4.11 d) 

k 

mvm 

s 

24i^ k Sin$ k 


(4.11 e) 

k 

mhm 

a 

24^ k 0os^ k - 24 SinjZf k 


(4.11 f) 

k^, 

BBT 

3E 

J2* k -0.5 Sin 2# k 


(4.11 g) 

M * 

rt 

Mm/R 
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k da = K + °* 5 sin 


*11 






M g 4# 

R 



i<^ (cosjz^ - cos#) d# 

R 



!L Cos# d# 
© 


M d# 

a 

E 


(4.11 h) 

(4.12 a) 

(4.12 b) 

(4.12 c) 

(4.12 d) 


I = — — (4.13 a) 

a t 

I? = J_ (-1—) 3 <4.13 b) 

a t 

a = width of the radial sectionof the arch 

t = depth of the radial section of the arch 

Eqns. 4*7, 4*8, and 4,9 have 3 un know ns ,V, H and M* and hence 
will give unique solution to determine the values of and 

putting those values on eqns. 4*4, 4,5 and 4*6 M, Rand S^. can 
be determined. 
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To derive the eqns. 4.8* 4.9, 4.10 the contribution of 
Sj to strain easrgy is neglected and the radial sectional properties 
of the arch are assumed to be constant with the variation of jZL 
The external loading is also taken symmetric about Z-axis • 

4.4 ARCH AffiM 

For the arch analysis for the proposed method 
the elementary strip of the circular arch is considered and this 
elementary arch strip is considered to be supported on the edge 
beam. This support is assumd to be elastic support which can 
defence- in three directions, i.e., vertical, horizontal and 
rotational directions „ In general each support point has three 
degrees of freed ome while actual conditions are determined from 
the boundary conditions of the -shell. The support reactions are 
determine d by satisfying the compatibility of the deformations of 
shell edge and edge beam at the interface in mid span. 

The elementary arch strip is analysed by the method 
described in article 4.3* The edge beam is assuad to offer 
restraint against vertical deformation only. Again considering 
the' statical equilibrium the arch is in equilibrium under the 
action of external loading and specific shear. 

As the arch strip in separated, for arch analysis, from 
edge beam then it would be in equilibrium under the loading 


of ]? y and V 
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Pig. 4.2 STATICAL EQUILIBRIUM OP ELEMENTARY ARCH STRIP 

n 

* « S’ 0 - statical equilibrium 2 p + V : = 0 

n 

or 7 « - Z P y (4.H) 

i =4 ' 

Therefore f the arch analysis essentially stands as the 
solution of tv70 oinultane ous eqns. 4.9 and 4.10 satisfyiig 
the conditions given by eon. 4 . 14 . Hence solution may be 
achieved solving two simultanious Equations. Putting the 
the value of V from eqns. 4 . 14 to 4.9 and 4. 10 

(a?3k mh + Tk rh } H + (j3j Srihrn> H ’ = E s\ * T \l 7 


(4.15) 
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W - - E s w m - T\ r (S 5 ^)? 
Where, 0? = B/t 


( 4 . 16 ) 


She unique solution of eqns. 4.15 and 4.16 will determined the 

values of d and h . As soca-as are known Mg and 13 g at 

any longitudinal section, are given by: 


M, 


% 


% + S V 7 +Rc nil H 


(4.17) 


M/ SS I + C V + f> F 
P s nv w c nh tl 


(4.18) 


where e 


ET 


Sin$ k - Sintf 


(4.19 a) 


! nh ~ 

Cosj^ - Qosjtf 

(4.19 b) 

! bv = 

- Sin$ 

(4.19 c) 

nh = 

- Cos$ 

<4.19 d) 


4.5 BOUNDARY COljDlflQIB 

Par single barrel shell H = 0 and M = 0 as it is ass ua s d 
that edge bean does not offer aiqy restraint towards horizontal and 
rotational deformations . 



lor the interior shell of. multibarrel group ^ = 0 ard 
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4*6 SUIVMY OF CQHHIFAfIGHS 

Far the system tic and easy computioml work of cylindrical 
shell analysis, the proposed method is detailed below in 6 steps* 


1. Sectional Iropertie s (R e£ .Pig .2 .3 ) 


a) Assume geometrical properties of the transverse section* 

b) h = R (1- CosJ^ k ) + d 

c) l c = 2R Sin^ k 

d) Area of cross section, k Q = 2 Rt# k 

A b = 2 bd 


e) Centre, of gravity. 


f) Moment of inertia.. 


A s = V *b 
! e =H 
Z b * h - a./2 

h = < A c Z c + A b Z b > / A a 
I„ = B 3 t k T 

C I 


\ = bd 3 /6 


Z rc= V Z c 


J rb‘ 


Zy~ 
b s 


h " V + V Z rt/ + *cKJ 

g) Co-ordinates of nodal points (£ef. Fig .2 . 4 ), 
dff = 0^/k 

k is always even number such that d0 > 5° 

lumber of nodal points n = m+3, m=k+1 

z ± =a(1-Cos^ i )-2 s 

Y^ = R Sinj^ i 
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h) Distance of top fiber from C.G., Z„ r+ = Z 

. yt s 

Distance of bottom fiber from C.G., = h-Z g 

i) First moment of area, i) For shell portion 


Qi = 2 Et (Z s - R k zi ) 
ii) For beam portion 


«i - b ( 2y b - 4 ) 


2. Loading (tfef. Fig. 1 -1c) 


Determine q u and 




4 


L. 


Sin/, 


(cig + t ) + q u ] l c + 2bdj 


3 . Beam Analysis 


a) T x t * ^ (2yt/ ) S'*" 1 * 

ff*b * «o ^ )2 (z yb /:i b ) Si ””’ X 

b) »•<>/* > <’* 

c) Specific shear = - ( q / 2 I ) /. 

V S 1 

4. The above computation canbe done systematically by the 
following tabular farm. 


1 ~ Co-ordinates 

Nodal — n r tir ~zy rr — qr 

point degree rad . 1 1 x 



Specific 

shear 


3.. Arch Analysis 
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a) Shell 

ds 4 - ds ffi * Bd//2 , ds t = Rd/ 
dy^ - R (Si^^- SinjZ^) 
dz^ = H (C os0^~ Cosj/^j ) 

where, V *V ^5~ 4 - ^ 

0j= 0±~i~ && i ~ 6,7 n-2 

4-1 = 4 = 0 


b) » Ti = q g ^ 

6x 


? hi = 


Beam 


is-j= ds 3 = d/4, ds 2 = d/2 


“ 0 

dz* = ds . 

1 i 


3=4,5 


VW + V “fV 


i i 

C j I Tj ) sin ^ - ( s 

3 = 4,5 3 = 4,5 


F hj ) Co aH t 


e) The arch analysis should he carried on with the help of above 
express ioniin the following tabular form. 


Sodal ds dy ds 
points 


dz l 0 d ^ \ ¥ l 


d) V 
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A. 


11 
A 12 

^21 

"^22 

*21 


R/t 


= *3 


k , + 5k, 
xnh nil 

tj\3 jj. 

Slam 


= A. 


12 

T'k. 


nm 


^ Iff 

12 / — s- 




R 


( Cos#, - Cos# ) d# 




*31 


4 

/ E_ Cos# d# 


’4 


k 


A 


k 


ML 


12 / JDSL d# 




R 


Hie above integral -ay be evaluated by 81mpson*s rule. 


Ch 


Or 


= E w- 


s 


*h - r *21 + **> 2 - ** W 


= E s " B - T 3 ^ - I 3 ^ 


IT 


A 


■] r n i rv> 


i ? 


! I as* ! ic. 


(4*20) 


H and M i Q solved putting bourdery conditions on 4.20. 


6 . Final values of 


% 


= 4- P 


P-C M J + Rc Ji + Hfir 


OJ 


ch 


- I\T + 0 7 + c„iH 

o nv rib 
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a e v £3 . es « kp V V,W W w ** 

arfc Slven ln appendix far different values of 0 and fL . 

J£ 

4.7 JCLLO ST5 A SPED SKAMPlg 


Example - 2 

Span = 30ru, 0 k = 300 , r « 9m, <1=9001;;, b =11.25 cm, t=9cm. 
Loading: Lime terracine a 45 kg/m 2 

Jjivc loc>d ss ytj kg/m 2 

Total =120 kg/m 2 

Step-1 

Sectional properties: 

h = 2.106 m, l c = 9.0 m, A g = 1 .0507 m 2 , Z q = Q.6466 m 
I s = 0.3795 m 4 

d^= 50 where k = 6, n = 7, n=lO 
z y t = 0.6466 m, Z yb = 1 .4594 m 

In shell portion Q. = 1 .62 (0.6466 -9k) 

In beam portion Q i = 0.1 125 (2.12 - Z 2 ) 

Step-2 

q ° ~ { x 3652.656 + 2 *0.9 x 0.1125* 2400)=4650.7061 kg/m 

Step-3 

Longitudinal stress at top fiber, - -717000 ljg/n 

longitudinal stress at bottom fiber, (f xb = 16.20000 kg/m 
0^0 = 58475*1435 

■' 4 * 

Specific shear = 6127.4165 Q ± 
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Step-4 

Table 4.1 STRESS RESUITAM1S ££, 


Toada. 

point 

t ' 

■ JZT 

Co-ordinates 

Q 

(Tx 2 

(kg/nr ) 

(kg/m) 

Specific 
shear 
(kg /m Z ) 

In In 

degree radian 

Z, 

(m) 

& 

(it?) 

1 



1 *4594 

4.5 

0 

1620000 

0 

0 

2 



1 .0094 

4.5 

0.1250 


7310.9125 

765.12 

3 



0.5594 

4.5 

0.2044 


11950.8231 

1251 .79 

4 

30 

0.5236 

0.5594 

4.5 

0.2044 


11950.1215 

1251 .79 

5 

25 

0.4363 

0.1966 

3 .8036 

0.2571. 


15020.3541 

1574,53 

6 

20 

0.3490 

-0.1 018 

3 .0782 

0.262 9 


15380.3618 

1610.05 

7 

15 

0.2618 

-0.3399 

2 .32 94 

0.2308 


13500.7451 

1413*46 

8 

10 

0.1745 

-0.5099 

1 .562 9 

0.1699 


j9?20:4175 

1040.50 

9 

5 - 

0.0872 

-0.6123 

0.7844 

0.0897 


525 0.3950 

549.95 

10 

0 

0 

-0.6 466 

0 

0 

-71 7000 

0 

0 


fa ®r? \ m mi&ymoiYO 2*> 
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1 0 

v= z F v; j = 356.0932 
0=4 

The shell is sirgle barrel, so,H = 0 and M* = 0 
Step - 6 

Pinal stress resultant is calculated and given in table 4*3 
and companed with the result calculated by classical theory . 

Table 4.3 AMD 


jef 

C rm7- — 

% 


— C y-. 

% 



mv 

1 

2 

IN 

1 

2 

30 

0 

0' 

0 

-0*5 

- 484 

- 180 

20 

0*15798 

121 

137 

-0.342 02 

-2583 

-2270 

10 

0.32635 

-110 

- 96 

-0*17365 

-4565 

-4257 

0 

0*5 

-424 

-460 

0 

-4934 

-5010 


Column 1 and 2 of and stands for the values calculated by 
proposed method and classical theory respectivty. The variation 
of stress resultants are shown in figure 4*3. ‘ 

4*8 DESIGMOOP REISFQRCEMEKES 

— — ~ ' 

a) Theoretical Aspects 

The 4 main stress resultants, for the purpose of design 






“600 -400 


-200 


too -5000 -4000 -jeoei -p«o iooo 6 


p^. 4.3 y/t RjATiON OF S~fg.£SS fcESUiTAWrj 
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of a cylindrical shell roof, are at x = 0, 1^,$^ and 3^ at 

x - i/2. On the basis of these stress resultants tie reinforcements 
are devided into two groups, 1) Force reinforcements^) Moment 
reinforcements » These are explained below. 


1 ) Forces Reinf orcements 

i) Longitudinal reinforcements, this is being to overcome 

ii) This is being due to 

iii) Diagonal tension created by the combination jof W x and 


ari( ^ a °tion on a small element of the shell, 

may be considered as orthogonal stress in a plane stress system. 



x -*3 1 ■ 

Fig. 4.3,, PRINCIPAL STRESSES III SHELL ELEMENTS 
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In true theoretical sense the ccmbinationshould be of 
where - ‘'l, 0 /p L - Bonnily is negligible 


Kfl = ^ 


B? 


_ X 

incipal stresses., ^ 1}2 = +* ^ (ET-: Ify) 2 +4 1^ j 2 . 


2 EL 


and the direction is given by tan 28. = 




V % 


(4.21 a) 
(4.21 b) 


Thus eqs . of 4.21 provide all the necessary information, to design 
the reinforcement, these should be examined at x = 0, l/4 and l/2 * 

At x =Ljk, 1^ = 0, therefore = H x and = ]^ 

where 9 = 0° or 90° 


At x 


= 1/4, H ! 2 = -1 

^ ri 


t % 


Vr 


end tan 29 =2 Sy / (2^- 1^) . 


± ^ at- b/-m 


X 


J 


At x = C> 5 1C. and IFv are zero. Thus S’. 


x 


■i,2 = i V 


and tan 2 9 =2 lT y ^/0 00 . 0 = 450 # 

• # 

2) Moment Reinforcements. This is transverse reinforcements 
to overcome the bendirg moment M^. 

The analysis of shell has been based on the assumption 
of homogeneity of the material. In other words, the elastic 
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20 m.m. 0 is provided fron bottom to shell 
upt 0 - 0 = 20°. 

ii) For Diagonal tension 

M a ximuia = 15380 kg/m. Area of steel required is 

10.9 cm 2 /m. Hovide 12 0 @ TO cm c/c . 
iii) For Transverse Force, 

is always compressive and maximum value is 5010 kg/m* 
So no reinforcement is required. 

2) Moment Reinforcements 

Moment reinforcement is required only f or®^f and it 
charges it’s Sign at 0 = 14°* 

TVEnrinmm = -424 Kgm/m. To resists this moment by a 

sirgly reinforced concrete section the eff. depth required/*? .49 cm 
if the eff. cover is kept 2 cm then eff. depth available is 7 cm 

* * 0 .K. 

♦ 0 

Steel required to resists this mount = 4-97 cm /m , provide 
1 0 / @ 1 5 cm c/c. 

Distribution steel: 10 / T @ 30 cm c/c. 

4.9 DISCUSSION 

a) The formulations, which have been derived in anticle 4*4 and 4*5 
and presented in 6 steps of calculations are fairly simple enough 
and can be eoneieve:&' without the knowledge of shell theory f with 
the help of the table s given, the computational work is simple and 
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less tedious. The tables which are given can be used for any 

values of R/t, l/l and 0^. The single barel, interior and 

eaterion shell of a mult ibarr el group can be analysed uniquely 

Stress resultant can be obtained at ary longitudinal section can 

^ith g.jn. n.TDc fry these formulations. The error involvement with the 
1 

calculations of arch analysis is less than ’that which are done' by 1 
column analogy method or elastic centre method. The error 
involvement ,with the problem of example ~3_ of article 4*7, for 
at top and bottom are 5.5% (more) and 6 .25% (less), for 
and errors are 15% (less), 8.5%(more) and 1.5% (less) 
respect ivly for maximum values 

b) The expression, of step 5(&) of the article 4.6 has been 
derived satisfying the static equilibrium. This could be done, 
by satisfying the vertical deformation compatibility. The vertical 
deformation of the edge beam at ary longitudinal section is given 


w. 


vb 


4 ■ % 

— -i-y ( 0.12319 V -0.19351 + JL 0.1231 9 bd f ) 


Sin 


n-jrx 


( 4.21 ) 


Erom the eqn. 4.8 the vertical deformation of the shell is given 
as (assuming external load varies with Sin mx) 


w. 


vs e 


~(T 3 k + Tk ) 2 7t)- (T 3 k • . ) H +(T 5 k )M +JpK,. 

v mv hv mvh mvm il 


i Sin ifcc 


( 4 * 22 ) 
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To satisfy the vertical deformation compatibility at the shell 
beam interface w vg should be equal to w yb . Therefore vertical 
deformation compatibility yields, 

(I 3 k n7 -KDk w -°.12319 Cfc) ^+(T 3 fc OTh )H+CT 3 k OTa )M'=1-a! 3 K 11 -tc b (w vbg -w vl:)s ) 

(4.23) 

Where, c fe = _ l4 , w ybo . = 0.12319 bdV , 

s bd 3 b 71 

d 


w. 


Vbs = °-19351 1^4 


The steps 5(d) of article 4.6 will be changed as follow: 
Step 5(d) of article 4.6: 

= 4 


t = R/t, <= b = fk _nl, 

bd- 5 


E„ 


w vbg = bd ,w yg =0 .19351 ^ 


A 11 = l 3 l w + Tk nr °- 12319 % 

A 12= A 21 = l3k mvli 


A 13~ ^31 ~ 


m 


^23 


2r k , + Tk , 
mh nh 


= A. 


A 33 - T 


32 

3 


^ ^rnh- 






1 


12 


M 


s Sin jZ^ d$ and , Kp 2 and are given 


\ E 


in step 5(d). of article 4.6 

Above integral may be evaluated by Simpson’s rule. 
- T Kji+ c b (w vb g- w vbs^ 

°2 - W ^*>1 + *22 


E -W - 
s m 


t3k 31 
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Therefore vertical, horizontal and rotational deformation 
compatibility at the interface yields, 



H 

M 



(4-24) 


The bouhdery conditions are: 

\ * 

i) For single barrel shell H=0 and M" = 0 

ii) For interior shell of multibarrel group, w^=0, w n =0 
substituting these boundery conditions m eq. 4.2 and solving 

I 

them 7, Hand M could be defined uniquly. The final stress 
resultants will be obtained from step 6 of article 4.^. 



CHART ER-V 


COIKLUSIOH AID REGOM'IEM)AT TOM' FOE FURTHER WORK 

The following conclusions are drawn based on the 
present investigation. 

5.1 EFFECT Oh SHELL DIMENS TOM'S ° 

1 . Reinforced concrete shells; Ifegnitudes of 
stress resultants decrease with increase of semicentral angle 
while increase with span to radius and depth to span ratios 
by both beam and classical theories. The influence of membrane 
solution on the total solution increases with semicentral 
angle. Beam theory agrees closely with classical theory 
for any values of radius to thickness and depth of edge .beam to 
span ratios when span to radius ratio is greater than 3. • 
Accuracy of beam theory depends on semicentral angle to a great 
extent. Higher value of semicentral angle for higher value 
of L/R ratio, yields the stress resultants with 4 to 10 % 
error. The preferable range of depth of edge beam to span 
ratio is 0.03 to 0 . 05 . 

Vertical, horizontal and rotational deformations 
increase with span to radius and depth of edge beam to span 
ratios while decrease with the increase of semicentral angle 
for classical theory. Beam theory also exhibits the same 
feature of vertical deformations. Restraints of horizontal 
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and rotational def or nations produce stress resultants which 
are in a nature to compensate each other. The ratio of the 
magnitude of vertical deformation, determined by classical 
theory to that determined by beam theory is in the range of 
0.262 to 1 .43 • Depth of edge beam to span does not have 
much influence as span to radius , ratio and semicentral angle 
have on this range. The range could be reduced as 0.8 to 
1 . 4 - 3 ,. with the proper selection of 0^. and Jj/r ratio. 

2. Erestressed Shell: Stress resultants are reduced 
by prestressing force for suitable depth of edge beam for 
all values of span to radius , ratio and semicentral angle 
by classical theory. Tension occurs at shell edge for depth 
to span ratio 0.05 and 0.04. Beam theory does not give ary 
tension for higher values of depth to span ratio. Ko exact 
conclusion could be drawn about the influence of semieentral 
angle and span to radius ratio on the agreement between 
classical and beam theories in the variation of stress resul- 
tants'. Kb shell with depth of edge beam to span ratio ^>0.04 
could be designed by beam theory. 

Beam theory underestimates the vertical deformation 

very much. From Table 5.1 it is evident in this respect two the 

a 

ories,. can not be compared except far/few cases. Horizontal 
and rotational deformation have higher magnitudes than those 
of reinforced concrete shell. 
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Table 5-1 COMPARISON OP BEAM AM) CIASSICAX THEORIES 


i/r 

4 

d/I 

e o 

e m 

' P 

e 

\ -■ 


6 

Reioarkj 

d 

d 

PSC’- 

ROC 

psc •' 

: RCC 


30 





31300 

5.5 

-0.135 

0.89 

3.33 

40 

0.03 

0.4 

0.4 

0.4 

53 ;50 27.75 

0.0715 

0.388 


45 





61 ,52 

45.82 

0.105 

0.262 


30 





30.00 

5.7 

0.078 

1 .030 

5.00 

40 

0.03 

0,4 

0.4 

0.4 

23-60 

10.0 

0.14 

0.800 


45 





15-00 

15.5 

0.16 

0.590 


30 





69.60 

8.0 

0.6 

1 .08 

10.00 

40 

0.03 

0.4 

0.4 

0.4 

69.2 

6 .2 

0.524 

1 .156 


45 





72.0 

11 -5 

0.458 

1 .43 



0.05*J 




56.5 

9.1 

0.016 

1 .510 

3.33 

30 

0 . 04* 

0.4 

0.4 

-3- 

* 

o 

50.0 

4 #0 

-0.123 

1 .250 



0.03 




31 .0 

5.55 

-0.135 

0.89 



0.02 




12 .2 

11 .5 

-1 .1 

0.73 




0.0 

0.0 


5-0 


— 





0.0 

0.4 


42.7 

- 


- 




0.4 

0.0 


4-5 

- 

- 

- 

3.33 

30 

0.03 

0.4 

0.4 

0.4 

31 .0 

- 

- 

- 




0.2 

0.4 


41 -0 

— 

- 

— 




0.4 

0-2*2 

40.0 

— 

— 

— 






0.8*4-25 -0 

— 

0.4 

— 






0.6 

’ 57.0 

- 

0.307 

- 

3.33 

30 

0.3 

0.4 

0.4 

0.4 

31 .0 

- 

-0.135 

- 






0.2 

14.0 

— 

-1 .77 

— 






0.0 


5.5 


0*89 


> = Maximum percentage of error in the variation of ^ by beam 

theory. 

> = j(w v ) beam theory/ (w v ) classical theory 

PSC and RCC indicate prestressed and reinforced concrete 
shell respectively. 

* 1 Tension occurs at the shell edge by classical theory 

< Vvfc U ■ 

2 At the bottom of^edge beam theory gives compression while 
classical theory gives tension 

3 Beam theory gives tension at crown. 
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5.2 EFFECT Off PRES1RBSS life ff ORGB : 

The reduction of stress resultants depends not only 
on the prestressing force but the combination of depth of the 
edge beam end and mid eccentricities. Difference between 
beam and classical theories decreases with sag. For shells, 
provided with axially prestressed edge beam, the beam theory 
agrees with classical theory within an error of 5 to 10 percent. 

Prestressing force is the most influencial parameter 
for the disagreement between beam and classical theories, 
w ith suitable combinations of depth of edge beam end and 
mid eccentricities. and decrease with the increase 
of prestressing force. Beam theory gives the same feature- 
In classical theory prestressirg has a very senative influence 
in the region of edge beam for the distribution of 11 b 

- influence in the distribution of , by classical theory, 

is not significant while beam theory indicates its.*. signi- 
ficant influence. ¥$ changes its* sign and decreases with 
the increase of prestressing force by both the theories. 

Shells with the coefficients of prestressing force ^0.4 
could not be designed by beam theory. 

The downward vertical deformation for non prestressing 
shell gradually decreases and ultimate attains upword 
deformation with the increase of prestressing force for both 
the theories. Rotational and horizontal deformations 
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increased with prestressing. 

Beam theory assumes loading is uniformly distributed 
over the surface or chord area of the shell and negleets the 
relative displacements within each transverse section. 

These are not correct for shells with edge beam. Bor 
reinforced concrete shell the corrective line loads for edge 
beam are not much while fear prestressed shell it has high 
values. The vertical displacements, estimated by beam 
theory for reinforfed concrete shell is well comparable 
with that by classical theory while for prestressed shell 
this is not true. Though beam theory assumes "that the rotational 
and horizontal deformations are negligible, which are not 
correct, but this assumption does not hare much influence 
on the agreement of beam and classical theories as the 
restraints against these deformations produces the stresses 
which are in a nature to compensate each other. The most 
influential point is the discrepancy in vertical deformation. 
This diLcripency increases with the corrective line loads 
due to edge beam, which inturn directly proportional to . 
the prestressirg force. This is the main reason of the 
disagreement between the shell and beam theory for reasonably 
higher values of prestressirg force. 
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5.3 PROPOSED METHOD; 

The proposed method which is detailed in 6 -steps, 

given article 4.6, has been formulated considering the 

variation of external load with Sin fzJL - in longitudinal 

1 

direction. The arch analysis portion, developed by strain energy 
method , is simple and less tedious. The errors in the 
computations with the help of tables which are given in 
appendix are small. Maximum number of simultaneous equations 
to be solved is 3 which is not'# all difficult for compution. 

She 13? with any type of boundary conditions can be analysed 
by this method. The values of the coefficients, given in 
the tables, are independent of structural geometry of the 
shell and. can be used for any values of span to radius 
ratio, depth of edge beam to span ratio, thickness to radious 
ratio and semicentral angle. Stress resultants at any longitu- 
dinal section can be determined by this method. 

5.4 SUMMARY OP THE COaCLUSlOIS; 

1 . Reinforced concrete shell with any values of R/t, 
d/L and can be analysed by beam theory with the proposed 
method for^/fe more than 3- 

Single barrel shells, exterior and interior of 
multibarrel shells could be analysed by beam theory with 
the proposed method. 
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•2. The error-' varies 4 to 10$ for 2 to 12 % 
for N^.0' Maximum error occurs for shells having higher 
values of L/R with lower values of 0^ and lower values of 
'0£ with higher values of L/R. 

Errors can he minimised with the proper selection 
of 0t£ as higher the value of L/R higher the value of 0^ • 

3. Erestressing has very sensitive influence on 
1^. in the edge beam while in the shell portion it is not 
3 i ch sensitive. Change of K ^0 with the introduction of 
pressing is not much significant. Influence of prestressing 
on 11 ^ and M 0 has same level of influence throughout the 
shell. 

4. The reduction of stress resultants is not 
dependent only on p hut on the suitable combination of 
d/L, e 0 /d and e m /d. along with p . 

Shells with axially prestressed edge beam can be 
analysed by the proposed method of beam theory. Shells 
provided with d/l^> 0.04 and coefficient of prestressing 
force more than or equal 0.4 can not be analysed by beam 
theory . 

6. The proposed method to analyse shell with beam 
theory is systematic, simple and less tedious. Tables, which 
are given to analyse a shell* can be used for any values of 
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R/L and R/t and 0 ^. The stress resultant at ary longitudinal 
section can be estimated by this method. 

5.5 SCORE OR FURTHER WORK ; 

1 . Analysis of 'cylindrical shells by beam theory 
with different support and leading condition. 

2 . It appears worthwhile to extend the analysis 
for the outer shells for 6 boundary conditions for beam 
theory and eight boundary conditions for classical theory 
and to study the effect of higher Eourier terms on the 
proposed method. 

3. Development of the algorithm of to attain the 
optimal design criterion of cylindrical shells by beam 
theory . 

4. It is worth-while to prepare the design table 
of prestressed concrete cylindrical shell by beam theory. 
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A. 1 SHELL COEEFICIENTS - 

Table A. 1.1 SHELL COEFFICIENTS 

r (0 Measured from' the Grown) 

F -2R (aB^-bBg) cos cosh a^0 - (aBg + bB^ )sin (3^0 sihh a^0 

(cB^-dB^) cos cosh a^0 - (cB^ + dB^)sin @J0 sinh a ^ 


F 

R 

B i 


B 2 


b a 

“jn 

- JOEL- sin kx 

X c 

n 2 (l+r) 


n 2 

n 2 ( V-1) 

n 2 


2Dk 2 sin kx 

1 


0 


0 

9x 

- 2T)k:^ nQ p kx 

1 


1 

-1 

1 

N* 

0 

■■J£E3f.t sin kx 

0 


1 

0 

-1 

. \ 

- 45rk 4 sin kx 

T 5 

.*•1 

1 

4- V" 

1 

1 - ^ 

9i 

- 2pk^ ROR kx 

V 

f + 2 


2 

V - 2 

2 

u x 

4Drk 5 eofl kx 

hEf 5 

-1 

1 

+ ^ 

1 

1 - f 

W 

r 

2 sin kx 

1 


0 

1 

0 


., with 

jf *=» 0, so 

that 

D * 

St 5 / 12 



k « me.. 
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Table A, 1.2 SHELL COEFFICIENTS (. ' from 'Oro 


A-2. 


^ 'SL 

T- *~ 

a * 

m m 

o o 

o o 

r- **t~ 

CCL CCL 

Jh ri 

*H *H 

CQ CO 


PQ CQ 
frj 

•fc + 

I" <? 

C6 O 


^ TS. 


*§ *§ 

*H *H 

CQ, CQ 

v*“ ***T" 
CCL CQ. 

CQ 02 

o o 

o o 


CM 

PQ PQ 

& n3 

T* lO 

R R 

<0 o 




A. 2 S O'BFFICI^TIS OF THB I'RAITS TBRS T! SRCTIOTO PROPERTIES 
Table A ,2.1 OOEPP'ICIEHDS FCR C .G- . 


0 

.in degree) 

. 0 

(in rad ) 

Cos 0 

Sin 0 

k z 

45.00 

0.78540 

0.70711 

0.70711 

0.09968 

42.50 

0.741 76 

0.73728 

0.67559 

0.08921 

40.00 

0.69813 

0.76604 

0.642 7 9 

0.07927 

37.50 

'0. 65450 

0.79335 

0.60876 

0.06988 

35.00 

0.61087 

0.81 915 

0.57358 

0.061 04 

32 .50 

0.56723 

0.84339 

0.53730 

0.05277 

30.00 

0. 52360 

0.86603 

0.50000 

0.04507 

27.50 

0.47997 

0.88701 

0.46175 

0.03795 

25 .00 

0.43633 

0.90631 

0.42262 

0.03143 

22 .50 

0.39270 

0.92388 

0.38268 

0.02550 

20.00 

0.34907 

0.93969 

0.34202 

0.02018 

17-50 

O. 3 O 543 

0.95372 

0.30071 

0.01548 

15.00 

0.26180 

0.96593 

0.25882 

0.01138 

12.50 

0.21817 

0.97630 

0.21644 

0.00791 

10.00 

0.17453 

0.98481. 

0.1 7365 

0.00507 

7-50 

0.13090 

0.991 44 

0.13053 

0.00285' 

5 - 00 

0.08727 

0.9961 9 

0.08716 

0.0012 7 

2 .50 

0.04363 

0.99905 

0.04362 

0.00032 

O 

O 

•* 

O 

3 

-0.00000 

1 .00000 

-0.00000 

-0.00000 



Table A.2.2 COEFFICIENTS FOR MOMENT OF IHERTIft. 
(0 Measured from crown ) 


0k 

k I 

0k 

k I 

* 

4 

k I 

K 

k I 

6 

45.00 

0.01216 

44.50 

0.01152 

44.00 

0.01091 

43.50 

0.01032 

43.00 

0.00976 

42.50 

0,00922 

42.00 

0.00871 

41.50 

0.00822 

41.00 

0.00775 

40.50 

0.00730 

40,00 

0,00687 

39,50 

0.00647 

39.00 

0.00608 

38.50 

0.00571 

38.00 

' 0.00536 

37.50 

0,00502 

37.00 

0.00470 

36.50 

0.00440 

36.00 

0.00411 

36.50 

0.00334 

35.00 

0.00358 

34.50 

0.00334 

34.00 

0.00311 

33.50 

0.00289 

33.00 

0.00269 

32.50 

0.00249 

32.00 

0.00231 

31.50 

0.00214 

31.00 

0. 001 98 

30.50 

0.00182 

30.00 

0.00168 

29.50 

0.00155 


* 


0k , given in degrees 



Table A. 3.1 COEFFICIENTS FOR ARCH CALCULATIONS 
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29.50 4.03178 0.11006 12.35693 - 0.52356 6.08484 - 1.06324 0.08629 0.94346 
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Table A. 4. t COEFFICIEN TS FOR STRESS RESULTANTS ( M } 



c vn 

°hn 

4 

c vn 

c lm 

45.00 

-0.70711 

-0.70711 

42.-50 

-0; 67559 

-0.73728 

40.00 

-0.64279 

-0.76604 

37.50 

-0.60876 

-0.79335 

35.00 

-0.57353 

-0.81915 

32.50 

-0.53730 

-0.84339 

O 

« 

o 

o 

-0.50000 

-0.86602 

27.50 

-0.46175 

-0.88701 

25.00 

-0.42262 

-0.90631 

22,50 

-0.38268 

-0.92388 

20.00 

- 0.34202 

-0.93969 

17.50 

-0.30071 

-0.95372 

15.00 

-0.25882 

-0.96593 

12.50 

-0.21644 

-0,97630 

10.00 

-0.17365 

-0.93481 

7.50 

-0.13053 

-0.99144 

5.00 

-0.08716 

-0.99619 

2.50 

-0.04362 

-0,99905 

0.00 

-0 ..ooooo 

-1 . ooooo 





Table , 

A. 4. 2 COEFFICIENTS 

; OF STRESS 

RESULTANTS (M^) 

‘ 0 

c vm 

c hm 

0 

c vm 

c hm 

0 k = 3Q-.00 



0 k= 32 . 5 O 



30.00 

o.. ooooo 

. 0.00000 

32.50 

0.00000 

. 0.00000 

25.00 

0.07738 

-0.04028 

30.00 

0.03730 

- 0.02263 

20.00 

0.15798 

- 0.07367 

25.00 

0.11468 

-0.06292 

15.00 

0.24118 

- 0.09990 

20.00 

0.19528 

-0. 096 30 

10.00 

0.32635 

- 0.11873 

15.00 

0.27848 

- 0.12253 

5.00 

0.41 284 

- 0.13017 

10.00 

0.36365 

- 0 . 14-142 

0.00 

0.50000 

- 0.13397 

' 5.00 

0.45140 

-0. 13-280 




0.00 

0.53730 

-0.15661 
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0 

c vm 

°hm 

0 

c vm Cfcm 

= 35.00 



$ c = 37.50 


35.00 

0.00000 

0,00000 

37.50 

0.00000 . 0.00000 

30.00 

0.73580 

- 0.04687 

35.00 

0.03518 - 0.02580 

25.00 

0.15096 

- 0.08716 

30.00 

0.10876 - 0.07267 

20.00 

0.23156 

- 0.12054 

25.00 

0.18614 - 0.11295 

15.00 

0.31476 

- 0.14677 

20.00 

0.26674 - 0.14634 

10.00 

0.39993 

- 0.16566 

15.00 

0.34994 - 0.17227 

5.00 

0.48642 

- 0.17704 

10.00 

0.43511 - 0.19145 

0.00 

0.57358 

- 0.13035 

5.00 

0.52161 - 0.20284 




0.00 

0.60876 - 0.20665 

= 40.00 



$£= 42,50 


40.00 

0.00000 

0.00000 

42.50 

0.0000 0.00000 

35.00 

0.06921 

- 0.05311 

40.00 

0.03280 - 0.02877 

30.00 

0.14279 

- 0.09998 

35.00 

0.10201 - 0 . 0818 ? 

25.00 

0.22017 

- 0.14026 

30.00 

0.17559 - 0.12875 

20.00 

0.30077 

- 0.17365 

25.00 

0.25297 - 0.16903 

15.00 

0.38397 

- 0.19988 

20.00 

0.33357 - 0.20242 

10.00 

0.46914 

- 0.21876 

15.00 

0.41677 - 0.22365 

5.00 

0.55563 

- 0.23015 

10.00 

0.50194 - 0.24753 

0.00 

0.64279 

- 0.23396 

5.00 

0.58843 - 0.25892 




0.00 

0,67559 - 0.26272 




$ c = 45.00 


45.00 

0.00000 

. 0.00000 

20.00 

0.36509 - 0.23259 

40.00 

0 . 06432 ' 

- 0.05894 

15.00 

0.44329 - 0.25382 

35.00 

0.13353 

- 0.11205 

10.00 

0.53346 - 0.27770 

30.00 

0.20711 

- 0.15892 

5.00 

0.61995 - 0.28909 

25.00 

0.28449 

- 0.19920 

0.00 

0.70711 - 0.29289 


